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^f~) • Various integrals over elliptic integrals are evaluated as couplings on spheres, resulting in some integral and series 

7—i ' representations for the mathematical constants n, G and ({3). 
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^ Introduction 

• ^ I Multiple elliptic integrals are integrals over expressions containing complete elliptic integrals. There are various 
/\f . techniques for computing multiple elliptic integrals, including the methods of Bessel moments 13,13, geometric 
transformations of Watson type 1 4], hypergeometric summations |5] and Legendre function relations (a). 

In this work, we showcase several geometric methods to evaluate multiple elliptic integrals in terms of the circle 
constant n := 4i:^o(-l)^(2/+ l)-\ Catalan's constant G := I^gC- 1)^(2^ + l)-2 (7,%, Apery's constant ((3) := L%-l£^ 
among others. We will make extensive use of techniques for harmonic analysis on spheres and hyperspheres llOll . 
in order to provide geometric interpretations for a large assortment of multiple elliptic integrals, new and old, famous 
or obscure. 

To flesh out, we will establish a series of identities by computing certain integrals on x (ID or x (31) 
in two ways: (1) applying (ultra)spherical harmonic expansions to the integrands; and (2) introducing (hyper)spherical 
coordinates with trigonometric parametrization. In 311 we follow up with some new perspectives on the Beltrami 
transformations introduced in our earlier work |6], in conjunction with applications of Beltrami transformations to 
selected multiple elliptic integrals derived in § ^1121 Some integral formulae in §3T]|3]will pave way for our subsequent 
works that frequently call for multiple elliptic integral representations of the constants n, G and ((3). We wrap up the 
main text with a collection of hypergeometric summation formulae (HJl that descend from certain integral identities in 
§ ^1131 In the appendix, we tabulate all the integrals and series proved in § ^1141 
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1 Integrations on x 



1.1 Spherical Harmonic Expansions and Geometric Parametrizations 

We will compute some integrals on the Cartesian products of two unit spheres x S^, using both the spherical 
harmonic expansions and geometric parametrizations. By evaluating the same integral using different methods, we 
are able to establish some interesting identities involving Catalan's constant G and the complete elliptic integral of the 



first kind K(^) := JqH - ¥■ sir? 6)-^''^ AO. 

Lemma 1.1 (Some Integrals on x S^) (a) Let c be the unit sphere equipped with the standard surface ele- 
ment da such that the surface area is m(S^) = fg2 da = 4n. For any unit vector n e S^, the following identities hold: 



8 f%(.2£ + l)^ 
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S2 4n |n-»i| |ni-n2l l»2-»l 
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8nVu(l-u)(l-u + uw)V'i--vwVw(l-v) 



vj n Jo 1-k 



j^iog-dk. (2) 



(b) Let P e {0,7i) be a given angle, and define the Heaviside theta function as 9h(x) = (1 + ■^^)/2,x £ IR \ {0}, then we have 
the following integral formulae: 
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(2/+l)/3 



f da,r 

Js2 471 Js 



i-'W):- ) ^cos 



do-2 



0ff(cOSjS-IIl ■II2) 



f'dq Cdrf' 

Jo Jo Jo 
n dq 
8 Jo v^l-gv^ 

r da,r 

JS2 471 Js 



s2 471 |» - ni I ^2(cos jS - m ■ »2) I'^a - »l 

d.- 



87iVa - r)(l - s)rs(l -q)Vq + tanHli/2) 



tan2()S/2) 
d(T2 1 0if(cosjS-iii ■112) 



= Cdq CdrC 
Jo Jo Jo 



S2 47t |» - ni I ^2(cos p-ni-n2)\f^2+n\ 

1 



ds- 



(3) 



C0S(jS/2) f^^l Ln . 2^] 

j Kh/A;2cos2^ + sm^ 



87r\/(l - r)(l - s}rsa - q) Vsec^ipm - s[qr + tan2(;S/2)] 

dfe=-icos?i r^^^s^dt 



cos(/3/2) 



Jrn/2 
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OdO 



2)Jo l + 2t^cosP + t'^ 



sinOV 1 - cos^ 6 smHp/2) 



(4) 



Here, the function 



1 rf^ 

£(/3) = G+- logti 
4 Jo 



7i-a , 
tan — : — da 



has special values £(0) - G, £(27i/3) - 2G/3, and £(27r/5) - £(47r/5) - 2G/5, leading to the identities 



G 
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2 Jo 8 Jo 
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dk 
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(V5 + l)^^K|^^^^f^\/A;2 + 5-2^Jd/fe-(^-l)^^K|^^^^^\//fe2 + 5 + 2^ 



dA; 
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Proof (a) Using spherical coordinates for points on the unit sphere n(9,(f>) - {sm9cos(p,sm6sm(p,cos6) e S^, and the 
famihar spherical harmonic expansion (with Pf being the Legendre polynomials, Y^m the spherical harmonics) 



1 ^P({n-n') ^ f Ye,nie,cl>)YemiO',(l>'} , , , „ ^ 



47i\n-n'\ ^0 4jr e^om=-e 2£ + l 

we may compute 



tU' js2 471 4jr In- ml |«i-«2l Iw2-«l hm=-e {2£ + lf to(2£ + l)^' 

r dCTi r dffa 1 1 1 ^ y^;„(0,0)y^;„(7r-0,0 + 7r) _ ^ (-1/ 

js2 4;r js2 4;r In-ml iMi-nal I"2 + «I hm=-e (2£ + l)^ to(2£ + lf' 

Here, we have exploited the orthonormal properties of the spherical harmonics 

sin0^d0j / d(l)jYei,n'(0j,cf>j)Yem(6j,(pj)^S(,/,5,n'm, J = l,2 

J-Ji 

to complete the integration with respect to »i £ and n2^S^, and have noted that 



4jr > -Pf(n-n), \/n,n eS , 

m^-e 2£ + l 

which equals 1 if we choose n' -n for Eq.[TJ and equals (-1)^ if we choose »' - -n for Eq.[2l 

In lieu of the spherical harmonic expansion, we may evaluate the surface integrals by variable transformations. We 
take, without loss of generality, the north pole n - (0,0,1) on the sphere, and employ spherical coordinates for m and 
»2, to compute 

r 1^ r dg2 1 1 1 

^ Js^ 471 Js2 An Im -#121 Iw2 -«l 

r ^i^fl rd0ir 02, „ rd02 1 

= / cos — dui I I cos — d02 / ■^^^^^^^^^^^^^^^^^^^^^^= 

Jo 2 J-n 'in Jo 2 J-jt An v^2[l - cos0i cos02 - sin0i sin02 cos(</)i - 02)1 

1 r du) r'' cos % cos ^ 

= - ^ d0i dflz , ' ' 

4 Jo n Jo Jo V2(l-cos0icos02-sin0isin02cos(p) 

If we make the substitutions 6i ^ n - 9i and n-62, the denominator of the integrand will remain unchanged, 
while the numerator becomes sin ^ sin ^ instead. By the addition theorem of cosines, we are led to the expression 



pd^ /"'d02^= 

^ Jo 8jrjo jo 



61-62 

cos o 



V2( 1 - cos 01 COS 02 - sin 01 sin 62 cos (p) 
Similarly, we have 

r dw r r sin 

Im^ \ d0i d0o ^ 

'-' Jo OTT Jo Jo V2(l-cos0icos02-sin0isin02cos(p) 

We then introduce the new variables (see Remark|(T)]below for geometric motivations) 



o(p 01 + 62 sin 01 sin 02 sin^ 

u = cos — , i; = sm , w-1- - 



2' 2 ' sin2^ sin2^' 

with Jacobian determinant 

Sin(psin|0i-02l 01 + 02 ^ / jz -; — ; r 

cot -2\/uw(l-u)(l-v)(l-vw). 



, ^ d{u,v,w) 
det- 



5(01,02,¥') 

If we swap 01 and 02, or perform an equatorial reflection on both polar angles n- 61,62^ n- 62, then the values 
of {u,v,w) will not change. It is then clear that the region < 0i < 7r,0 < 02 < 7r,0 < < tt is mapped to the unit box 
< u < 1,0 <!;<1,0<m;<1 exactly four folds and we have the geometric correspondence 

01-02 . 01 + 02 r- 

cos — - — = vl-i;if, sin — - — -sjv, 1 - cos0icos02 - sin0i sin02COS(p = 2!;(1 - u + uu;). 
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All this allows us to put down 

Jm= [ du [ dv [ dw — 

' Jo Jo Jo 87i\/u{l-u){l-u + uw)\/vw{l-v) 



8jr\/u(l -u)(l-u + uw)Vl- vw\/w(l - v) 
Now, we note that the complete elliptic integral of the first kind satisfies 

•'O v^l-/fe2sin2 Jo 2x/u(l-u)(l-k^u) 
After integrating in v, we may reduce the triple integral representations into expressions involving K(\/l-f2)- 

i Cd. Cdi r ^ . i Cd, [\ f ^Kfyr^lde; 

T-i 4Jo Jo Jo ji^ua-u)a-u + u^^Wva-v) 4Jo Jo ^u{l-u){l-u + ue) 2Jo ^ J 

421= i f'd. f'd^ f ' ^-^ . J- r 

^ 4Jo Jo Jo jiVua-u)a-u + u^^Wl-v(^\/T^ ^ttJo 



1 f^K(./T^) 1 + f 

■log- — -d^. 



Here, the complete elliptic integral in question satisfies an identity called Landen's transformatiorQ [Ref 11 , p. 39] 
which is equivalent to the following observation of GauB: 



I(a,b):- j — = 7^1 . ,Vnh 



Va^cos^e + fo^sii?^ ^ 2 
Therefore, we may simplify with a Mbbius transformation t]-{1- 0/(1 + 0- 

4j= 7— dry, 42^=-/ log- 

^ Jo 1 + 77 nJo l-T] 7] 

which completes the proof of all the claimed identities in part (a). 



(b) To prove the series representations for Eqs.[3]and|4l we need to avail ourselves with the Mehler-Dirichlet formula 
(see item 8.927 in UJ] or §4.5.4 in flil): 

~ e Yem(ei,,pi)Ye,n(e2,<p2)cos[(e +\)fi-\ ~ 0H(cos;S-ni -na) , , , , ^ 

47r2^ 2^ -— ; — ^2^Pe{ni-n2)cos[(e+^)fi]^—= '^ :, |«il = |»2l = 1- (5) 

e=am=-e Zff + 1 ^=0 V2(cos/3-ni -112) 

To justify the triple integral representations, we introduce the change of variables (see Remark [(1)] below for geo- 
metric backgrounds) 

cos/3-cos(0i + 02) cos/3-cos0iCOs02-sin0isin02COS(p l-cos(0i-02) 
1^ ^ — ; ' '■^ 2 77^ — TT-^ ' s = 



cos/3+1 cos/3-cos(0i + 62) 1-COS01COS02 -sin0isin02cos^ 

with Jacobian determinant 

2sin0isin02 sin(0i + 02)sin(0i - 62)sm(p 



d(q,r,s) 

det- 



3^1,02, cp) 



(C0S/3+ l)[cos/3- cos(0i + 02)](1- COS01COS02 - sin0isin02 coscp) 
2 v/(l - - s)s(l - q) Vsec^(/3/2) - s[qr + tan^(/3/2")] + tan^(/3/2) 



^The simplest way to verify Landen's transformation 

' ' •^0 v/l-(l-f^)sin20 •'0 v/co?eTf2~n2i J- 

tanfl + sec 



v/co?eTf2~n2i V(l + f)^-(l-f)2sin2fl 1 + f vl + f 

is to use the trigonometric substitution 

6 = arctan 



4 



As before, we may verify that the image of the region 0<6i<n,0<92<7t,0<(p<n covers the unit box <q < l,0<r< 
1,0 < s < 1 exactly four times, and there are geometric relations 



61-62 



sec^ — -s\qr + tan^ — I cos — , sin 



cos- 



coSjS - cos 01 COS 02 - sinSi sin 62 cosip - 2qrcos —. 



q + t&rr — cos — , 
2 2 



Consequently, we infer that 



71(71- (3) d(p 0if(coS|S-cos0icos02 -sin0isin02COS(p)cos 



v^2(cos P - cos 01 cos 02 - sin 01 sin 02 costp) 
1 



r£^rd0ird02 

Jo 071 Jo Jo 

Jdq I dr I ds — 
■'0 Jo 871 V(l - r){l - s)rs{l -q)Vq + tan2(/3/2) 
71 dq 

^^0 y/l-qyyq + tan^(j3/2)' 



81-82 



1 rP 7i-a 

I[4]:-G+—l logtan da- j dq l dr j ds — 

4 Jo 4 Jo Jo Jo 87t\/il 

Jo Jo 



4:7iVr(l-r)(l-q) 



r)(l - s)rs(l - q) VsecHpm - s[qr + tanHp/2)] 



cos^W gr + tan^ ^ ) . 



In the last Une, we may replace r hy = p :- qr e (0, q), to obtain 



cos(/3/2) /■! 



471 



f'dpf , ' K 

Jo Jp Vpiq-pXl-q) 



COS — A/p + tan^ — 



d9: 



cos(/3/2) /■! 



KK//fe2cos2^ + sin2^^ 



d^. 



Now, with the integral representation for the complete elliptic integral of the first kind, we have 

cos(/3/2) . 2 



- f d^ f dk 

Jo Jo V(i7?y 



2 - uHk^ cos2(jg/2) + sin^(i3/2)] 
fi] {l + t^)logt 



1 

- / arcsm — — 
2jo ^ 



2icos(/3/2) df 



cos 



2jjo 1 



dt. 



. + 2t^cosl5+t^ ^ 

after integration by parts in t. The last line of Eq. |4]is evident once we can verify the integral formula: 

-71/2 

0<k< l,O<0<7r/2. 



K V/^^ (,Qg2 _|_ gjjj2 



^ 



-k^cos^eV l-sin2 0sin2 
To show this, we pick an arbitrary non-negative integer n, and adapt the foregoing proof as 



r KiJk^cos'^(P + sm^cl)]k''dk^ C dq f\r f^ds — 

Jo I ; Jo Jo Jo 47i\/{i- 



iqr) 



n/2 



r)(l- s)rs{l - q) \/l - s{qr cos^ + sin^ (p) 



(6) 



Such a transformation is valid so long as the integrals converge on both sides of the equation above. Then, we take a 
"box-to-box transformation" (cf Remark [(5)l below) 



1-qr 

with Jacobian determinant 



which results in 



det 



d(q,r,s) 



1-W 



f^K{Jk^cos^(P + sin^ (p]k''dk^ C d'^ C dY C dW ^ 

J<^ \ I Jo Jo Jo 471 V^yrd rWl - y{W cos^ + sln^ (p) 

dr dW — = / d0 1 dk—— 

Jo 4 Vrrd -y)\Ji-yw -y(i-w) sin^ ■'o Jo 



- k"^ cos^ sin^ sin^ 



5 



Here, we have set T - sin 0/(1-^ cos 0),!^ - k in the last step. As we run through all the non-negative integer 
powers « = 0,1,2,... , we see that Eq.[6]must hold. 

The special values taken by the function 2,{fi) result from the following integral relations for Catalan's constant G 

nnli /■7I/12 /■3h/20 

-2G^2 Iogtan0d0 = 3 / Iogtan0d0 = 5| logtanSdS, 

Jo Jo J;i/20 

which in turn, are consequences of factorizations for tan 30 and tan 50 into products of tangents. ■ 
Remark (1) In the proof of the lemma above, we introduced the following change of variables 

^(p _. 2^1 + ^2 _ sin0isin02 _ sin^^i^ 



U - cos — , V - sin , W -1- „ n n - ^ n n 

2' 2 sin2^ sin2^ 

cos/3-cos(0i + 02) COS /3- COS 01 cos 02 -sin 01 sin 02 COS (p l-cos(0i-02) 
9= 3 — ; ' 7, 77. 7r. ' 



cos/3+1 cos/3-cos(0i + 02) 1 -COS 01 COS 02 -sin 01 sin 02 COS 

as a wholesale argument. Now we illustrate the line of thoughts that underlay these geometric transformations. 
Taking the integral 

■ 61-62 



\/2[l-cos0i cos 02 - sin 01 sin 02 cos (p] 
as an example, we may naturally reparametrize with u - cos^ f , which brings us the expression 



61-62 J 
COS 2 d» 



nyjuil- M)(sin2 - u sin0i sin02) 

We now introduce the new variables 



01 + 02 

V = sin , W = sin0isin02, 



with Jacobian determinant 



det 




det| 


3(01,02) 







2 COS 2 2 2 

^cos 01 sin 02 sin 01 cos 02 



01 + 02 

cos — - — sin(0i - 02) 



= Vi - y2\/y2 _ ■jyy'i + - y2. 



Such a map sends the square region < 0i < 7r,0 < 02 < tt to the parabolic triangle 0<y<l,0<W<y^ exactly four 
folds, because the values of y and W will remain intact if we trade the point (0i,02) for (02, 0i) or replace the point 
(01,02) with {71-61,71-62). This allows us to further transform Jjjjinto 



JfYI=4x — f^dy r dW f^ — 
^ 16 Jo Jo Jo n\ 



du 



VT^y/V'^ - W^ud - ii)(y2 - uW) 

- f dv f dw f — 

8 Jo Jo Jo 7rv^u(l -u){l-u + ui 



iw)Vvw(l - v) 

where we have set u = y^,^; = (V^ - W)fV^ in the last step. 

If we apply the foregoing arguments to part (b) of Lemma ll.ll we may deduce 



7i(7i-p) 1 r'' r cos^ig^du 



Id Jo Jo J(0,l)ni 
8 JJJo<u,v,w<l;vil-u + uw)>- 



dudvdw 



r~r-, x.cosfi-i 



2 7i\Ju(l-u)( — I ^-v-uv + uvw)\/w(l - v) 

due to the geometric constraint ni ■»2 < cos/3. We next move on to further reduce the triple integrals, starting with two 
observations: (i) An inequality j;^sin2^>sin2 S follows from the geometric constraint 

cos /3 > cos 01 cos 02 + sin 01 sin 02 cos ^ > cos(0i + 02 ); 
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(ii) The introduction of a new variable 



t;u;-sin | _ cos/3 - cos(0i - 62) 

o '.— a — ~ ^ 

cos/3-cos(0i + 62) 



-,1 



with 



det 



d(S,u,v) 



diu.v.w) 



,2/3 



leads to a factorization 



^ \-v-uv + uvw = [1 - u(l - S)] iv - sin^ ^ 



It is thus reasonable to introduce another new variable r :- 1- u(l - S) e (0,1), so that S < r is a consequence of 
u = (l-r)/(l-S)< 1. 



With the replacements 



u = 



— , y^gcos — hsin — , 

1-S ^2 2 



1-S oP + tan^ £ 

W-S+ sin —- —, 

V 2 ^ + tan2§ 



namely, a transformation of variables 
V - sin r 



cos^ § 

with Jacobian determinant 



£(0,1), r=l-[i(l-s): 



uud-u;) 
1 -^e(0,l). 



,^„2/i 



tan^f 



-,r 



diq,r,S) 



d(u,v,w) 



vHl-w) [q + tan^]{l-S) 



we may convert the triple integral in li,!;,^; into 
^ 7i{n-p) 



[v-sin^) 
) 



gcos^ ^ 



^r(l - r)(r -S)a-q)Jq + tan^ f v/ gS + tan^ § 



Finally, it would be natural to rescale S as 



S + 



tan^ 



r+ ■ 



tan^i 



1 - U + UW 



so that s ranges from to 1. 

(2) The "sphere-to-box transformations" (0i,02,<p)" {u,v,w) and {9i,62,ip)^ (q,r,s) used for proving Jjjj Jjg} ^ ^iid ^ 
are extensions of the techniques employed in Watson's work on simplifying certain triple integrals of physical 
significance. The spherical harmonic expansion trick is a continuation of the author's own previous research |17] on 
the Hilbert-Schmidt norms of certain operator polynomials in the mathematical model of electromagnetic scattering, 
which involved an integral on x x x S^. 

During the course of our proof, all the interchanges of integrations can be justified without theoretical difficulty: 
the integrands are non-negative and the triple integrals are convergent. 

(3) The major purpose of the lemma above is to build confidence in a geometric method, rather than producing results 
that are particularly surprising. Although we have brought a variety of formulae under the framework of xS^, some 
of these results can be derived via alternative methods independent of spherical geometry, which in turn, can be traced 
to certain standard references. 



The integral identity 



1 K(r?) ^ jt^ 

+ 7] ' 8 



is well known (cf item 615.09 in 1 11] or item 6.144 in 1 12]) and can be proved by simpler means invoking only double 
integrals. Our integral representation for Catalan's constant 



log 



d^ = G 



7 



is equivalent to formula 4.21.3.13 in which can be derived using identities satisfied by hypergeometric func- 
tions 119], without going through the integral on spheres. 

The identity 

7i{7t-p) _n r'^ Aq 

8 sJo y/l-qs/q + tanHpi2) 

is elementary, while the equality 

wa\-G+-\ logtan da = - cos— / ^ jAt 

4jo ^ 4 \ 2)Jo l + 2t^cosp + t^ 

can be double-checked by differentiating both sides in /3. An equivalent form of Eq. [6l as given by I. Mannari and C. 
Kawabata in 1964: 

■ KWC + S-CS), 



f 

Jo 



JO \/(l-Ccos2 0)(l-Ssin^0) 
was mentioned in Setting jS = in 



^ cos(/3/2) r"'^ BdB 1 71- a ^ 

G = — ItL^ I / logtan da, 0<p<n, (7) 

2 Jo sin0Vl-cos20sin2(/3/2) 4 Jo 4 

we obtain 



1 S i-n, 

2 Jo sm0 

which is a familiar integral representation of Catalan's constant (see 3.747.2 in 1,12.] ). 

(4) The variable substitutions in this lemma, either used in isolation, or in conjunction with hyperspherical geometry 
(see Proposition [2^, can facilitate the integration of many other expressions involving the complete elliptic integral of 
the first kind K: such integrals can be thus represented by certain mathematical constants, or converted into integrals 
of elementary functions. As a stand-alone example, we point out that the set of transformations employed for proving 
Eq.|4]can also be used to justify the following identity: 



V 2(cos p - COS a) J t^'^ '^^d [ [ ^ 0ff(cos/3-»i ■ra2) 1 

J/3 2jicosHp/2) 4 js^~^Js^ 4:71 \n-ni\ cosHp/2) |n2 + w| 

- [ dq f dr f ds ^ 

Jo Jo Jo 4;rv^(l - r)(l - s)s(l - g) V^l - s[qrcosHp/2) + sin2(jg/2)] 

^ j^K^]jk^cos^^ + sin'^^ kdk, 0<p<7i. 

In particular, for P-0, the integral on x obviously evaluates to unity, thereby bringing us the result 

1 r" f a\ 71 — a 

-I sin- log cot da= / K(^)^d^ = l. (8) 

jrJo ^ 2^ 4 Jo 

(5) We may write Catalan's constant G as a triple inte^al using either the specific case of /3 = or Jjg} Namely, 
the number 47r^G, as a member in the ring of periods 12011 . has two integral representations 



Jfl pi pi roo 

dq dr ds dt 
Jo Jo Jo 

J^oo rl rl rl 

dt du dv d 
Jo Jo Jo 



^2 + (1 - r)(l - s)rs{l - q){l - qrs) 



1 

w- 



+ u(l - u){l -u + uw)il - vw)wil - v) 

We can pass from one quadruple integral to the other by the following algebraic mappings between two unit boxes: 

w 1-r 

q-v, r-l-u + uw, s- u- , v-q, w-rs, (9) 

1-u + uw 1-rs 

fully oblivious to the geometric interpretations of the variables u,v,w and q,r,s. 

More generally, evaluations of many multiple elliptic integrals amount to the confirmation of relations among cer- 
tain members in the ring of periods, which is conceivably accessible by algebraic means |20]. It is sometimes possible 
to transition from one multiple elliptic integral to another using a variety of (geometrically motivated) birational map- 
pings. The Beltrami rotations and Ramanujan rotations in [Ref Hi §3] are examples of such transformations. □ 
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Upon integration, Eqs.[3]and|4]give rise to some integral formulae involving Apery's constant ((3) - ^ and 

complete elliptic integral of the second kind E(^) - f^^^il - sin^ 6)^''^ AO: 



8 16 



dcTi C dcr2 



JS2 471 Js 



S2 4:71 In -nil 



-K 



/l-rai-n2 



n l+\/l-q Aq 



|W2 



- — -— f log — 
-n| 8 Jo 1 



^/T^^-q 



" 7r(:fr -3) 6 
cos 



4 Js2 471 Js2 4jr In-; 

~ 6 Jo 



g / 2V(l-rarra;)72 
l+V(l-nrn2)/2_ 



»il V^l+V(T^ni-n2)/2l'»2-nl 



3-P'2 



1 1 3 

5 V 
4' 4 



1-9 



dq 



cos - d p = - : 

8 2 ^2 



dai f da2 1 n 1 



r doi r 

JS2 471 Js 



S2 471 In -nil 2 |n2 ■ 



[K(vT^)-E(vT^)] 



(l_^)3/2' 



dfTi /■ dcr2 



Jo 4 Js2 47r Js 



S2 4?! In -nil 
ds 



-K 



/l-ni-n2 



Jo Jo Jo 

= / log:; — ilogf 
Jo 



2 
1 



|n2 + n| 



47rsv^(l - r)(l - sMl - g)(l - gr) 



arctan 



(l-q'r)s 
1-s 



1 + i 



dcTi r dcr2 



/■;i/2 

Jo 



sin0 COS0 



■d0, 



r£(/3)cos-d/3 = 2v^log2= f — f 
Jo 4 Js2 47r Js 



g / 2V(l-rarra^ 
' ^' l+V(l-nrn2)/2 j 



S2 47r In-nil v^i + _ . „2)/2 l«2 + «l 
j'^^'M (f+i)tan-^^^ ^-(i-j)tanh-^ 



nil - sWil-r)rsil-q) 



3F2 



/ 1 1 3 

2'-^' 2 


(l-gr)s] 


5 7 
V 4' 4 


1-s J 



t-i 

dai f da2 1 7r 1 



t+i 



logtdt 



[ —I 

Js2 An Js 



S2 47r |n-nil 2 |n2 + w| 



.1 M .1 (l-g.s)E(/S)_(l_,)K(/<| 

j dg / dr / ds- 

Jo Jo Jo 

-ir^*^^^^^=/o 



qr)s 



qrs j 



4n(l - qr)s\/ (1 - r)(l - s)rs(l - g)(l - grs) 

nl2 



E(cos0) 
[ sin0 



sin0K(cose) 



QdB 
cos^S 



In the above, most integrations are straightforward, except a few that additionally require Taylor expansions. For 
example: 



" cos(/3/4)d/3 
Vg + tan2(/3/2) 



cos(/3/4)cos(/3/2)d/3 



COi 
Jo y/l 



-(l-g)cos^(/3/2) 



v/27r^(-ir(l-g)"- 



= £(-l)"(l-g)" 

n=0 ■'O 

(2m+1)! 4v^ 



„ /"^ v^cos(/3/4)[cos(/3/2)]2"+M/3 



n\Y{\-n) 



n=o «!r(i-7i)r(«+|) 3 



3-F2 



/" i 1 3 
2'-^' 2 

5 1 
V 4' 4 



1-g 



In [Ref.ll, Lemma 1.2], we employed the Hobson coupling formula to establish the following integral identities for 
01,02 e[0,7r): 



1 r^'' I 0\ fK sin^ K sin% , 0i + 02<?r 

- Ksin-d(/)=^ 

4-^0 I 2j [k cos^ K cos| , 0i + 02>7r 



4 Jo 



2sin(0/2) ] 



l + sin(0/2)^ 



d0 



v'l + sin(0/2) 



1 -^1' / 2sin(ei/2) 

Vl+sin(0i/2)Vl+sin{02/2) IV l+sin(0i/2)J 



^l+cos(ei/2)Vl+cos(02/2)^(v l+cos(ei/2) 



2cos(ei/2) 



2sin(e2/2) 
l+sin(02/2) 

'\^{ I 2cos(e2/2) 
j^l^V l+cos(02/2) 



(10) 



(11) 



, ei + 62>n 



9 



where cos - cos 0i cos 02 + sin 0i sin 02 cos (/>. In the light of this, we may use Eqs. \T0\ and [Tl] to "decouple" the integrals 
on X in Eqs.HDIlEIlllDandllE!] With some elementary trigonometric identities that were used in the symmetric 
reduction of ijj] and ijg] at the initial stages of the proof for Lemma [Lll we arrive at some integrals on two-dimensional 
simplices: 

- = -ff Kfsin^lKfsin^lcos^i^d0id02, dSQ) 

16 2jr JJo<ei<7r,o<02<Ji,o<ei+02<Ji v 2j \ 2j 2 



j^l I 2sin(ei/2) Urf / 2sin(e2/2) 1 

27r JJo<ei<7r,o<02<Ji,o<ei+02<7i \/l + sin(0i/2)\/l + sin(02/2) 2 

— — = — K sm— K sm— sm — — d0id02, ®]) 

4 JJo<ei<:ff,o<02<7i,o<ei+02<7i V ^/ V ^ 



(HQ) 



/ 2sm(0i/2) Urf / 2sm(e2/2) 1 
_ „ 1 rr ^ V l+sin(ei/2) ^ V l+«n(02/2) . + ^2 , . , „ 

2V21og2=— — —sin — - — d0id02- 

27r JJo<ei<7r,o<02<7r,o<ei+02<jr \/l + sin(0i/2)\/l + sin(02/2) 2 

The integration of the relation 

o.«^ n 1 . COS()8/2) /■! / 7^6 TT^K ^ 

£(/3) = G + - I logtan da= ^- — / K\\ k^coB^-+s\n'^- Ak 

4jo ^ 4 2 h W 2 2] 

was omitted from our hst of identities in Eas. [3ni3P^l4F1 and l4r**[ but will be treated in W^\ 

1.2 Integrals for Azimuthally Anisotropic Couplings 

In the physical parlance, we have computed in 31.11 some Coulomb-like interactions on spheres where the "force 
potential" is modulated with polar anisotropy. While reducing the triple integrals in part (a) of Lemma ll.ll we inte- 
grated over u before w; we then introduced an additional factor in the integrand (which modified the isotropic Coulomb 
potential) dependent on »i •»2 to the triple integrals in part (b). 

In the next proposition, we will redo part (a) of Lemma ll.ll bv completing the integration over w before u, to derive 
identities of a different flavor, and will explore the effect of azimuthal dependence of the "force potential modulation". 

Proposition 1.2 (More Integrals on x -S^) (a) We may rewrite Eqs. Uland^for any unit vector n e S^: 

Aai r do'2 111 



^ - f f 
8 js2 in Js 



S2 471 |wi-ra2l Iw2-»l 

1 



G 



Jo Jo Jo 87rv'w(l -u)(l-u + uw)\/vw{l - 1;) 
_ 1 r'^ logTjde _ logging de _ r'^ logtanfde _ ri logtdt 

4 Jo COS0 Jo COS0 Jo 2cos0 Jo t^-l' 

_ r dai r da2 1 1 1 

JS2 4:71 Js' 



S2 471 |wi-ra2l Iw2 + »l 

1 



871 V u(l-u)(l-u + uw)\/'i-- vwVwil-v) 



= f\u f\v Cd 
Jo Jo Jo 

-^r '^''11' -2 i\ikw{VT^dkU^ login d,. 

27rlJo v^TT^ Jo ^ I ] 271 Jo ^ ^1-^ ' 



v^l^ 

In particular, Eq.^^implies the identity 

VK.(k)fd k 



C'^^-2\\ik)^^dk. 
Jo yiTp^ Jo V ) 



(12) 



(b) More generally, suppose that the function a)(cot^) is continuously differentiable in (pfor 0<q)< (po, where q)o e [0,7i], 



10 



then we have 



a){k/Vl-k^)dk 
.(<Po/2) ^{l-k^)[l-k^a-i^)] 



1-^ e 



w(cot^^i^)0ff((po-|0i-</'2l) 1 



|W1 -»2l 



|ra2 + «l 



271 Jo [J cos 

r dCTi r dg2 1 — 2 

js2 471 js2 471 Ire -rail 
1 [K{k)fa){k/Vl~^)dk 

2jI jcos(<po/2) \/l-/fe2 

1 p r dcjiK)] f\r\. 1 , Vi-^^t^-^VT^ 

/ (d(k) + k— djc / dt I d^ , , log ; — 

TT Jcot(|<pol/2) [ dJC J Jo Jo ^tVl + K^t^Vl + K^(^t^ ^"f 

^ 2 ; 2 Jo Jo (t\/l + cot2(|(po|/2)\/l + <2^2 cotHWo\m 



log ; , (13) 



where 



rai-ra2-(rai-ra)(ra2-ra) . 

101 - 02h= arccos ^ e [0,7r] 



V^l-(rai-ra)Vl-("2-n)2 
Some particular cases ofEq. \13\lead to the identities 



f 

Jo 



Jo 



TT^ 1 dk 



[K(^)r- — ^ — = — log2 
[K(^)]2 



7C(3) 



4 ) ^ 2 "° 4 ' 
^ 7r ri[K(jfe)]2^ 



, log — dk - — \ — dk. 

VT^ Wi^ 2 Jo 



(14) 
(15) 



Proof (a) Integrating over w first, we may deduce 

1 



Cdu Cdv r 

Jo Jo Jo 



dw- 



Pi pi logi^du 2a 1 r 

■t;) Jo JO 87rM\/(l-u)i;(l-i') 4 Jo 



d77 



87r\/u(l - u)(l - u + MK;)-\/t;it;(l 
Meanwhile, it is not hard to verify 

8 ^ Jo 1 + 7? Jo [Jo a + rjWi-rj^^in'^g 
thus Eq.[l0follows from the simple algebra 

^2 '-/2iogl±£2||d0 /--/nogd + cose) 



1 ^;i/2 1ogi±«>s|d0 
/ = l-cos6 



COS0 



d9 



rJl/2 

Jo 



log(l + cos0) 

COS0 



dB, 



^ _ 1 f 
"8"~ 2 Jo 



COS0 



Jo 



COS0 



de 



-L 



'^''^logsinede 

COS0 



and 



8 Jo 



71/2 log .ff^de 



rnl2 

Jo 



logtanfde /■llog^d^ 



_ plog^d^ 
~Jo l^TT' 



Before treating the triple integral for Eq.[2B we make the substitution u - k^/{1 + k^) and w-t^, so that 

di; 



~ 271 Jo VTT^Jo [Jo Vi + K^t'^VT^ 



We then move on to manipulate the bracketed integral in t. It is straightforward to check that the bivariate function 

dt 



Jo 



satisfies the following partial differential equation 



f{K,v) Kdf{K,v) df(K,v) 
■ + — ; h U- 



2 dK 



dv 



(16) 



11 



with "boundary conditions" 



Urn v[f(K,v)f^O, lim f(K,v)^— __, , 

lim K[f{K, v)f = 0, lim xlfix, v)f^O for < y < 1. 

K — 0+ K — +00 



for Q<K < +(X>; 



Therefore, we obtain after integration by parts 



f(K,v) K df{K,v) df{K,v) 
+ — ; 1- V 



2 dx 



1 

G=- dx dvf(x,v) „ . „ , . ^ , 
71 Jo Jo [ 2 2 ax ov 

.±r\K{-^]r^-irdxf\u\fJ^]'--^r'^^ 

2nJo [ \s/TTt^]\ 1 + k2 TiJo Jo [ 2 \ 2% Jo 47rjo Jo 



Here, the last double integral can be further simplified, as we write 

2 



Vf{x,v)Y^ 
and integrate in x and v separately: 



di 2 r^^^ t 

Jo ^/T~^Vi + x^tA ~ J° J° ^/T^Vi + x^t^ Vi - vih^ Vi+K^et^ " 



C \r[nx,v)fdix\Av^4.l^ At Ca^ 
Jo (Jo J Jo Jo 



lOS" 



Then we may integrate by parts with respect to t to deduce 



1^'"' w 



tH ^ i-f 



1 1 + ^ 

KCO-Trrlog 



2f 



Therefore, we can verify Eq. [20 with the computation 

— dx dv[f{x,v)f^- K(OK 
4jr Jo Jo TT Jo 



1-^2 df 



■I 



1K(VW2), 1 + ^,- 

log di. 

2n^ ^1-^ 



Now that the last integral in the equation above evaluates to G (according to Eq. |2), we reach the identity claimed in 
Eq.m 

(b) We note that the transformation cos^((p/2) -u- x^/(l + x^) entails x - cot(|(/)|/2), and 

dffl r dCT2 1 W(coti^^i^)0if((po-|01-02l) 1 



r da,r 

Js^ 471 Js 



s2 471 Ira -nil 



= — f\v Cdt f 

2n Jo Jo Jet 



|ni -»2l 
a)(x)dx 



|ra2 + n| 



't(IV>ol/2) y/a + X^)(l + X^t^Wl - vfVT^ 



can be reduced to a double integral by an integration over v, and a subsequent transformation jc - k/V 
2ji Jo Ja 



ci){x)dx 



271 Jo [Jcot(lv,ol/2) y^{l + X^){l + X^t^) 



log 



1 + tdt 1 



1-t t 2ji 



Jo Jcos(ipo/2) 



w(^/v^^^)d^ 



V^(l-^2)[i_^2(i_^2)] 



log 



1 + tdt 



1-t t 



This establishes the first line of Eq. [13] Using integration by parts as in the proof of part (a), we may arrive at the 
expression 



1 [K(^)]2w(^/v^i^)d^ 1 r°° 

iTT Jcos(«)n/2) \/T^^k^ 47r Jcot(|ipol/2) 

2 



27r Jcos((po/2) Vl-^2 



(JL)(X)+X 



dw{x) 1 
d?c 



Ijc / dv[f(x,v) 
Jo 



w cot cot / / cot ,v\\ dv, 

471 [ 2 j 2 Jo [' I 2 ' j\ 



12 



where the last term vanishes for cpo - n, as in part (a). This is related to final identity claimed in Ea. 1131 as we note 
that 



Jo Jo Jo Jo ^/i 



t 



Jo Jo 



— „ „ ^ „ „ „ log- 



For 0<(po<7i, setting (i){k) - l/?c,cot((po/2) <k< +oo in Eq.[13l we obtain 



mk)fdk 1 + VsmHcpo/2) + cos2((po/2) 1 + ^df 
log — ^^ — ; — log- 



mk)]^dk _ ri 

Jcos(¥.o/2) k Jo"" (l + OcOS{(po/2} ( 

+ f'dtf' 

Jo Jo 



d( , =log- 



' ^t\/l + f cot2(|^o|/2)^l + f2i2 cot2(|(po|/2) 
which is equivalent to 



C i[K(^)]2_^ll^. ^^^^ 1 + V^sin^(.,o/2) + e cos^(yo/2) 1 + ^ d< 

Jcos{(/)o/2) 1 4: I k Jo 1 + f 

+ f dJ' 

Jo Jo 



„ „ , „ „ „ =log- 



' itVl + t^ COt2(|(po|/2) V^l + ^2^2 cot2(|(po|/2) 

since we have, by Taylor expansion and the monotone convergence theorem 



Jo i Jo to^n+l ( to(2n + 



,2«+l ( „to(2re+l)^ 
In the limit of (po — ► tt, we obtain 



As we may compute (using Taylor series expansion of the integrand, when necessary) 



Jo 



, 2 ^ 1 + fdC 7C(3) 
log log — log2 ; 



CdtrdAo^^^-'^-Cdtrdi'-io,^^-'^''^'' 

Jo Jo l-f Jo Jo ( 



1-i dt 

, r\ 2t\ogt /■i2arcsini 

= - dt df ^ ^- logtdt 

Jo Jo Jo v^T^ 

/■I (arcsini)2 tt^ 7^(3) 

= / dt- — log2 , 

Jo t 4^8' 

we arrive at Eq.[l4]as claimed. 

With a logarithmic input, we apply a similar limit procedure to Eq.[13]and write down 



fi^iog-^d.^r [ 

Jo Vl-/fe2 v/l-fc2 Jo Jo 



logJcdK 



Vl-^2 Jo [Jo v^(l + Jf2)(l + K2f2) 



, 1 + id^ 

log 



+ 2f (l + logK)djc/" d^ f d^ ^ , log 

Jo Jo Jo .,,/TT7272,/TTr2T2^ 



1 , v^rW-^vT^ 



Here, the integral involving logK can be evaluated in terms of complete elliptic integrals (see item 800.05 in or 
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4.242.1 in fl3l ). which brings us 

Jo VT^ 2 Jo V I ^ t Jo [ 2^ ^1 



1 + e 



K Vl-e^ 



^ ,K(v^^)iog(f^2) yrr^-^^iTTJ 



Jo Jo 



-f 

Jo 



1, i+f 



/o it^ 
|Ka)(2-logO-7log 



1-e 



Jo Jo 



K[Vl-^'|d^ 
log- 



The last double integral can be simplified via integration by parts 



1 /■! K(v/i^)iogi vi-efi-(VT^ 



Jo Jo 



d^ 



■log 



Jo Jo 



dfK \/l-f 



1 + logi 
t 



and an integration in t (using item 800.01 in or 4.242.4 in iflill for the term involving logO 
'i/l + logi 



Thus, we have 



L 



-K|Vl-^'| + Ka)log^ 



1 [K(^)]2 ^ J J. 

log , - nk 



=-['1 

2 Jo 



df 



2 Jo 







d^, 



as stated in Eq.[T5l ■ 
Remark (1) Using Tricomi's formula [Ref I21], p. 103] for the Fourier expansion of K(sin0), one may explicitly compute 



7T S^i [Tin+h]^ 



Jo v^iTl? Jo y I 4 to (nl)^ 



1111 

2' 2' 2' 2 
1,1,1 



as recorded in Eq. 73 of 11] and Eq. 35 of Igl- (Here, 4F3 is the generalized hypergeometric series.) This backs up our 
Eq. [T2I which was derived geometrically. 

(2) After a preliminary literature search, we have not been able to locate previous reports of the closed-form evaluation 
given in Eq.[14l though certain similar-looking integrals (Eqs. 76, 81 and 82 in had been considered in the frame- 
work of Tricomi's expansion, and the expression 7f(3)-27r^log2 arose in quite a few different contexts |22, 23, 24]. Later 
in this work, we will give some other integral expressions that evaluate to the same number as Eq.[l4l Furthermore, a 
"lower-degree analog" of Eq.[l4]that involves Catalan's constant is well known: (cf [Refill!, item 615.05], [Ref [13, item 
6.142] or [Refill item 2.16.2.3]) 



n 



71] dk 

K(^) - - — = -2G + 7rlog2. 
2 J k 



(17) 



(3) We may recast Eq. [15] into 



f 

Jo 



log- 



dk 



Jo 



[K{k)Vdk^—jF6 



I 1111115 
2'2'2'2'2'2'4 



drawing on a recent result of Borwein et al. mentioned in la], which asserts that 



2^\K{k)fdk^^^\K\y\ 



1-K^II dK^—lFo 



I 1111115 

2'2'2'2'2'2'4 

1,1,1,1,1,1 



with being the generahzed hypergeometric series. As pointed out in , the leftmost equality in the formula above 
is a direct consequence of Landen's transformation. □ 
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In 33.11 we will give a self-contained proof of Eq. [17] For the moment, we will accept the truth of Eq.[T7l and combine 
it with Eq.[l4]into another integral identity. 



Corollary 1.3 (An Integral Involving n, G and C(3)) We have the following formula: 



L 



k 



l--E(^) 

71 



7^(3) 

d^ - — 2G+ . 

4 2n 



Proof First, by virtue of the ordinary differential equation 



_d 

di 



til-t) 



dK(v^) 
dt 



K(Vt) 



(18) 



(19) 



we may apply integration by parts to the following integral: 

dK(v^)l 



[K{Vi)flogtdt^4 KiVDlogt— t{l-t) — -^\dt^-4 
Jo Jo dt [ dt I Jo 



d [K(v^)logf] dK(v^) ^ 
t(l- 1) — — — dt 



dt 



dt 



i-o+ I di J Jo di I dt \ 

= 71^ + 4 f^^^[E{Vt)-K(\/t)]dt+ (^[K(^/t)f\ogtdt, 
Jo t Jo 



which leads to 



Jo 



K(^) 



[E(^)-K(^)]d^ 



= -f' 

2 Jo 



[E(v^)-K(v^)]di=- — . 



(20) 



k ^ ^ 2 Jo t 

Then, as we add Eq.[20]to the left-hand side of Eq.[T8j before invoking Eqs.[T7]and[T4l subsequently we may compute 



Jo 



K(^) 



l--E(^) 

71 



dk- 



f 

Jo 
4 Jo 



1--E{k) 

71 



dk + 



71 2 r'^Kik) 



4 71 Jo 



[E(^)-K(^)]d^ 



l--K(^) d^: 

71 



71 r^K(^)+f 



4 Jo 



k 



l--K(^) 

71 



d^-2G + 7rlog2 



= --2G + 7^1og2- 
4 



2 

71 Jo 



,^,U)f-^-Y-^^^-2G.'-^, 
4 k 4 271 



as claimed. 



The reduction of a logarithmic factor in the integrand of Eq. [15] is not purely accidental. It can be viewed as a con- 
sequence of Parseval-type identities for the Tricomi transform (see FRef \26L §4.3, Eqs. 2 and 4] or FRef EtI. Eqs. 11.237 
and 11.241]): 



c){^f)ix)dx^O, 



f(x)g(x)dx - j\^f)(x)(^g)Mdx ^If^ \f{x)i^g}(x) + gix}(^f)(x)\ log^ dx, 



(21) 
(22) 



where f e LP{-l,l),p > 1; g e L'^(-l,l),g > 1 and ^ + ^ < 1, where the Tricomi transform is defined as a Cauchy 
principal value 



(^f)i 



x) := ^ / - 

J-l 71' 



fiOd^ 
ix-O' 



a.e. xe(-l,l). 



and induces a bounded linear operator ^ : L^(-l, 1) — ► L^(-l, 1) for 1 < p < +oo FRef lid p. 188]. In the next corollary, 
we explain this "Tricomi pairing" by putting Eq.[T5]in a slightly more general context. 

Corollary 1.4 (Some Applications of Tricomi Pairing) The following chain of identities hold: 

[K(v^l-jf2j] ^iT<-f^ [K{k)f 'ik^^j^ K(OK [ yr^] log ^ d(^^j^K{k)K[ Vl-k^] log ^ d^ 
1 r'^[K(k)f-[K{y/T^)f. 



-I 

71 Jo 



JO 

where arccosx £ F0,7r] for x £ F-1, 1]. 



■log^d^ 



jr Jo 



K(fe) K(\/l k ) - 2^^)d^, 



(23) 
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Proof In [RefM Corollary 3.2(b)], we showed that 

■1 K{Vl-K2)dK k ^ K(Vl-K2)dK 



K(^) 



j-i n(k~K) n j-i 



Q<k<l, 



which entails a Tricomi transform formula: 

■1 K(v^l-?c2)d)c k 



r^KW 

J-1 71 



ik-K) 



K(k), ^e(-l,0)u(0,l). 



(24) 



(25) 



Now, setting f(x) = gix) = K(v^l-x2) in Eq.lH and reading o&i^Oix) = (.^g)(x) = |fiK(x) from Eq.|25l one obtains 
the first equality of Eg. [23] 

Using the variable substitution ^ - (l-k)/{l + k) along with Landen's transformations for 2K((1 - ^)/(l + ^)) = (1 + 
k)KWl-k'^) and K(2\/k/(l + k)) -{1 + k)'K.{k), we may convert the rightmost term in the equation above into 

^ K(OK (yr^j log ^ de = ^ K(^)k( v^^) log ^ d^, 

thereby completing the first line in Ea.l23l 

To move on to the next term of Eq.|23l we fall back on the following formula (see [Ref d, Eq. 26] or [Re£ 17, Eq. 41*]) 



[K(v^)]2 = - r 

TT Jo 



2 /■! K(^)K(0^)dM 



0<t<l 



(26) 



to compute 



[K(^)]2-[K(vT^)]2 



-f 

71 JO 

7t^ Jo Jo 1 - 



log^d^ 



^1"^ ^iK(^) K(0^)dM 
/isin^e 



71 Jo 



[K(^)]2 ^ ^ ^ 

log - rife 



71 Jo 



logtanede 



K( ^/T^) log( 1 - M) d 

Vl-M 



[K(sin0)]^logtan0d0 



27r jo 



where the integration over £ [0,7r/2] can be rendered in closed form. A substitution k - ^ in the last step will 
bring the equation above to the last term in the first line of Eq. |23] In particular, this result also hearkens back to 
Eq.[15l 

Now, we note that the following Tricomi transform formula (see [Ref'^, Corollary 3.4]) 



/>( 



1-^ 



2d( 

7l{x - 



K 



1 + x 



K 



1-x 



-1 < X < 1 



allows us to specialize Eg. 121] into 



pj(x) 



K 



1 + x 



K 



1-x 



>dx= - 2 p{^f){x)K 



1 + x 



K 



1-x 



dx 



(27) 



for arbitrary / e L^(-l,l),p > 1. From [Ref 'is', Eqs. 12A.5 and 12A.12], we know that Vl-x^f(x) = logi^^ entails 
-arccosx, so Eg. [27] further specializes to 



1 [K(^)]^-[K(v^l^)]2 



71 Jo 



I 



\ogkdk ■■ 



1 arccosx^ 

271 J-l \/T^ 



1 + x 



K 



1-x 



4jr J-l 
1 

dx = — / 

71 Jo 



^1^ 



K 



1 + x 



K 



1-x 



>dx 



, arccos(l-2/5 )dk, 



which completes the second line of Eq. 
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2 Integrations on x 



2.1 Ultraspherical Harmonic Expansions and Geometric Parametrizations 

Now, we shall extend the geometric approach to integrations on the Cartesian product of hyperspheres S" x S", 
where n>2. The usual spherical coordinates have a natural extension to the unit hypersphere S" :- {x-{xi,... ,Xn+i) e 

' Xl- COS 01 

X2 - sin 01 COS 02 
X3 - sin 01 sin 02 cos 03 

Xn - sin0i---sin0„_icos0„ 
Xn+i - sin0i ■ • ■ sin0„_i sin0„ 



(0 < 0^ < 7r,j = 1,2, . . . - 1;0 < 0„ < 2n) 



so that the standard surface element (induced Lebesgue measure) can be put into 

do-= nsin"-^'0jd0^-. 

In place of the spherical harmonics Y/m, we need the Gegenbauer polynomials C^}\0 defined by the generating function 



= ^i^C</'(0, 0<i<l,-l<^<l. 



Lemma 2.1 (Convolution on Hyperspheres) Let S" c (n > 2) be the unit hypersphere equipped with the stan- 
dard surface element da and surface area m(S'^) :- fgnda - 2ji^'^'''-^^'^fr(^^^), then the product of the following three 
functions defined for n,ni,n2,n' £ S" 

^ 71-1 j ^ 71-1 j ^ 77-1 j 

, . ^ /■^^^Z ("-"i) , , ^SfCf' (ni-n2) ^ ^/C/ (n2-n') 

f(n-ni)^ fo+ 2^ , g{ni-n2)^go+2^ , hin2-n')^ ho+ 2^ 



will integrate to 



e=i n-1 ^ti n-1 ^ti n-1 



("-!) 

r dcri r do-2 , C"/ ^ (n-n') 

/ ,<^nJ ^"' ■ • n2)hin2 ■ » ) = /'o^o/io + L Z'/^^'i^ 

Js'^ miS"-) Js'^ niiS"-) 7^^ 



e 

(«-l)(2^ + «-l)2^ 



so long as both sides converge. Furthermore, the above integral formula remains valid for n-2. 

Proof To complete the integration with respect to »i e S" (re > 2), we exploit the addition theorem of Gegenbauer 
polynomials (see 8.934.3 in UJ] or §5.3.2 in 113]) in the form of 

(-1) ^ (71-3)! f 22A^-7)![r(^+j)]2 y2 72 

^'^^■"2)=^f^i3i^E in-2 + eljy. ('^-2 + 2j)[l-(ni-«)2]^ [l-(n-n2)Y >< 

><df^^n,.n)&,f'^n.n2)&'^'' ni-«(ni-n) «2-«(«2-«) 



'^~J '^~J J \\n-i-n{ni-n)\ \n2-n{n2-n 

We assert that only the j - term in the finite sum on the right-hand side of the equation above will eventually 
contribute to the integral over rai e S". This is because we may rotate wi about the axis n, without altering the value 
of » ■ » 1 , amounting to a latitude-longitude breakdown of the integral on hypersphere in the following manner: 

f(n-ni)[l-(ni-n) \ C " (ni-«)C. " ■ -■■ ^ — 7^ 

J J \\n\-n{ni-n)\ \n2-n(n2-n)\i m{S"-) 



L 

J-l ■' [J|s| = l,s-re=0 ■' I |«2 -»(«2 •»)! ; 



0, j=l,2,...,^. 



m(S«) 
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Here in the last line, dcTg stands for the surface element of S" ^, which forms the "equator" with respect to the "north 
pole" n; the integral on dag vanishes due to the orthogonality relation of Gegenbauer polynomials (see 8.904 in lll2n ): 

r"c*'^'(cos(p)sin"-2(/)d^= C d^\cos(p)&r^\cos(p)sin''-^(pd(p^O, n>2, j ^1,2,...,£. 
Jo Jo 

Now, to compute the surface integral 

dcTi „, _ ^ r da 



Js" m(S") Jsi m(S"-) 



~ {n-2)W. geCf^ (wi-n)C^" (ww2) 



n-1 



we employ the orthogonality relation for the Gegenbauer polynomials (see 8.904 in 111211 ) 

nin-2 + £y. 



2«-2/!(^ + ^)[r(i^)]2 



Especially, with the latitude-longitude decomposition, we have the relation 



r dai 

Js" m{S"-] 



which proves the identity 



Cy '{n-ni)Cy '(ni-n)=^^^^^|^C*/'W/''(0(l-^^)^de = 



miS" 



(n-2 + e)\ n-1 
e\(n-2)\ 2£ + n-l 



>e'e, 



r dCTi 5^ ' 

/ ,n,„J (n-ni)g{ni-n2)^fogo+ l^ffg^T 
Js" m(S") in 



C\ ^ (n-n2) 



l)(2^ + n-l)' 



so long as both sides exist in the distributional sense. Yet another application of the addition theorem for Gegenbauer 
polynomials will lead us to the formula 



f f ■^7^Jin-ni)g{ni-n2)h{n2-n')^fogoho+Y,fegehe — 

Js"- rnKS"-) Js"^ miS"-) (n 



C) ^ \n-n') 
-l)(2^ + «-l)2^ 



as claimed. It is worth noting that for « = 2, we have precisely 



C/ (n-n')^Pe(n-n')^4n ^ 



Y(,niB,cp)Yemiff,cl,^) 
2£+l 



for n{6,(p) - (sin0cos(^,sin0sin0,cos0) £ -S^, n'{6' ,<p') - (sin0'cos(/)',sin0'sin0',cos0') e S^, so the current lemma is 
indeed a generalization for the technique of spherical harmonic expansion employed in Lemma [Lll ■ 

Remark The integral formula in this lemma remains true in the "n 1""" limit". After we recall that, for every positive 
integer £, the Gegenbauer poljmomial has limit behavior 

(— ) 

C. ^ (cos0) cos£(h 
lim , 

71-1 £ 

we may immediately formulate an analogous integration formula on xS^: the product of the three functions defined 
for n - {cos (p, sin (p),n I - (cos0i,sin0i),»2 = (cos(^2,sin02),«' = (cos(/)',sin0') £ 

, ^ ^ ^ cos^(0-</)i) , ^ ^ cos^((/.i-02) , ~ COS^((/)2-0') 
f(n-ni)^fo+l^f( -p ' g^ni-n2)^ go+t.ge z , h{n2-n)^ho+ l^hg 



integrates to 



/ -T— \ ——fin-ni)g{n\-n2)h{n2-n)^fogoho+l^feg(h. 
Js^ 2n Js^ 2n 



cos£{<p-(p') 



Obviously, the truth of such a statement can be established by the elementary addition theorem of cosines and the 
orthogonality relation for Fourier series, which are just degenerate cases of the addition theorem and orthogonality 
relation for Gegenbauer polynomials. □ 
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In the next proposition, we narrow down to the case of x S^, where the unit hypersphere is customarily 
parametrized by n(tfr,6,(p) - (sini//sin0cos(^,sini//sin0sin(/),sini//cos0,cosi/'), along with the standard surface element 
do- - (sin^i//di//)(sin0d0d0) and total area fgsda - /o'^sin^i/zdi// /g'^sinede f^" A(P - In^. For n-2>, the Gegenbauer 



polynomial C^f'iO is identical to the Tchebychev polynomial of the second kind U^iO, satisfying 



C^"(cos0)=C7^(cos0)=<( 



£+1, = 

sin(^ + 1)0 



, O<0<n; 
sm0 

(-1)V+1), 0^71. 



Proposition 2.2 (Some Integrals on x S^) (a) Let c be the standard unit hypersphere, on which two unit 
vectors n e and n± £ are orthogonal to each other (i.e. n ■ n± - 0), then we have an integral representation of 
Catalan's constant 

r dai r dao 1 1 1 If l + sinw 

Js3 2^:'' Js3 271"' |n-nir |»i-n2r l»2-raj_r 4:Ji Jo l-smi// 

Generally, if two unit vectors neS^ and ne £ are separated by an angle (i.e. n-ne- cos0) where < < n, we have 
the identity 



lm[U2{ke'^)] _ ^ k^~'^ sm£0 _ r dai r dCT2 1 1 1 

1 \ ksinw If, l-cos(w + 0)l , , , ^ 

arctan — log diff, 0<k<l, (29) 

sm0Jo ^ 



^sin0 

sini// 



arctan - 



27rfesin0Jo I l-^cosi// 
where 

oo 

Li2(z):=^^, \z\<l 

n = l 

defines the dilogarithm function for complex valued z eC. 
(b) We have the following integral formulae: 

^1 zK(VT^) 



l-cos(w + 0) 

log 



1 - cos(i// - 0) 



Jo 



r— 

Jo Ja 



V(l-z^f + 4z^t 



df = Li2(z)-Li2(-0), |0|<1 (30) 



a-z*mvr~t) i i+z 

:df=-log- , |2;|<1 (31) 



2 Jo 



[(1 - 22)2 + 4zh]y/{l - ^2)2 + 4^2f z \- Z 



K(^l^)log — dt = zLi2 (z) - 0Li2 (-Z) - (1 + z)log(l + 0) - (1 - 0)log(l - z), |z| < 1 



(32) 



where the univalent branches of the square root and logarithm functions are chosen conventionally. As a consequence, 
Catalan's constant can be represented as 



G=---log2-- K(v^)log^^ ^ -dt, (33) 

4 2 4Jo t 

1 /■'^^2 p/2 sin0(l-cos*0)sin(^cos(/)K(sin0)K(sin(/)) 

— I I 2 o o — Q« ' (34) 

^ Jo Jo (sin + 4cos2 0cos2(/))3'2 



n/2 />7t/2 ^r>ofln _ r.r>a4 



and dO j dcj) 



(sin^ 9 + 4 cos2 6 cos^ cf))^'^ 
Proof (a) Specializing Lemma [2TT1 to the case n-3 with 

fo^go^ho^l; fe^ge^he^2,\/£eZn [1, +oo), 



cos0(l-cos 0)sin0cos(/)K(sin(/)) 

(35) 
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and using the special value of the Gegenbauer polynomial C^^'(O) - 11^(0) - sin^^^j^, we may readily confirm the left 
halfofEq.m 



r doi r 

Js3 2n^ Js 



do-2 



4 sin 



(-1)^ 



is3 27r2 |n-ni|2 |„i-„2|2 |»2-«±|2 /bo(2^ + 2)2 ta'^2e+l)'^ 
where n ■ wj_ = 0. More generally, if we start from g( - 2k^ instead of g'^ = 2 for positive integers (, we may derive 

do-2 1 1 1 



Js3 2jr2 jga 



2^ A^-^sin^© 
£l ^2 sin© ' 



0<k<l, 



js3 2n:2 |ra-ni|2 |w]^-^W2l^ l«2-«0l^ 
for « ■ We = cos 0, which is the first line of Eq. l29l 

Alternatively, we can evaluate the integral on x with an admixture of ultraspherical harmonic expansion and 
geometric parametrization. Following the proof of Lemma l2.11 we compute 



do-i 1 1 _^ k 

's3 27r2 |n-ni|2 Im -fen2|2 fr'^ 



e-i 



sin£'F 



'sin>F 



^sinV' 



■ arctan 



1-^cos?'' 



where n n2 - cosf . We then complete the integration over W2 using the coordinates n - (0,0,0, 1), w© - (O,O,sin0,cos0) 
and n2 - (sini//2sin02COS02,sini//2sin02sin(/)2,sini//2COS02,cosi//2)- This amounts to a verification of the second line in 
Eq.im 



Im[Li2(fee'")] 
^sin0 



do-2 



■ arctan ■ 



^ sin 1/^2 



l-^cosi//'2 

2n 



fs^ 2n^ |ra2 -M0p^ sin 1/^2 

= — T / sin^i//2di//2 / sin02d02 / d02 
4?!^ Jo Jo Jo 

= — / sin2i/A2di//2 / sin02d02 
271 Jo Jo 



( 1 - cos 1//2 COS - sin 1//2 cos 62 sin 0)k sin 1//2 



arctan 



k sin7//2 

1-A cos 1/^2 



( 1 - COS tfr2 cos - sin y/2 cos 62 sin 0)^ sin y/2 



2nk 



1 

sin0 Jo 



arctan 



ksiny/ 
l-^cosi// 



log 



1 - cosd/' + 0) 
1 - cos(i// - 0) 



dy/. 



The particular case for ^ = 1,0 = 7i/2 also brings us the integral representation of G, as claimed in Eq.[ 
sinf - sin(7r - y/), we may reduce the integral formula into 



Noting that 



G 



1 f 1 + sinw 1 f'^ l + sinw 

— (7r-i//)log- : di//=- log:; : dy/ - 

In Jo l-smw 4 Jo 



47rJo l-sini// 4 Jo l-sini// 

which is equivalent to the statement £,{n) - for Eq. H) 
(b) We first verify Eq.[30]for z-ik where < ^ < 1, namely. 



1 r"'^^ n-2w ^ 1 71- a ^ 
-/ logtan — - — di//=--/ logtan— — da, 

2 Jo 4 4 Jo 4 



Im[Li2(j^)]: 



1 ikKiy/T^) 



:d(, 0<k<l. 



V(l + ^2)2_4^2^2 

To show this, we consider the specific scenario of Eq.[29]with < ^ < 1 and - 7t/2, which amounts to 



Im[Li2(jA)]: 



1 

2^ Jo 



arctan ■ 



ksiny/ 



l + smw , 1 r 

log dy/- — / 

1-smy/ 271 Jo 



J!/2 



l-^COSI/' 

We now integrate by parts in y/, with the knowledge of Eq.|4]in Lemma [Lll 



arctan 



2k sini// 



l + sinw , 

log- : dy/. 

1 - smy/ 



Im[Li2(j^)]: 



2 f''^ 
71 Jo 1- 



Jo 



-2^2(.(,g2i// + ^4 

Then, we set ( - V1-k^ cosi//, to convert the integral above into 

fcosv ^K(^l^) 



cosi// 1 K \/l-(l-jc2)cos2i// dK 



dy/. 



2 /■'^'^ 

Im[Li2(jA)]= - / - „,2 o 74 



2 /■'''2 ^(l-^2)j,Qg^ 
7t Jo 

2 

TTJO 1- 



-2kHl-2rf) + k'^ 



I/; 



\/ C0s2 1// - ^2 
'^ CK(v/l-f2) 

V^l-r?2-^2 



de 



di// 



df 



drj : 



1 f^KCv^l^) 



Jo ,/(T 



v/(l + /fe2)2_4^2^2 
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where we have completed the integration in ry e (0, to arrive at the last expression. 

We then note that both sides of Eq. [30] are complex-analytic functions of z in the simply-connected domain |z| < 1, 
because the roots of the quadratic equation (1 - Z)^ + AZt - are given hy Z+- l-2t + 2i \/t(l-t), satisfying |Z+ 1 - 1. 
By the principle of analytic continuation, we may extend the identity from the top-half vertical diameter {z - ik where 
< < 1) to the open unit disk (l^l < 1). With the Vitali convergence theorem, we may then check that both sides of 
Eq.[30]are continuous up to the boundary \z\ - 1. In particular, the special cases z-1 and z = i for Eq.[30lrecover a pair 
of identities we showed previously in Lemma [Lll 



2jG = Li2(j)-Li2(-j)= / -c\f^i I K(^)d^ = 2iifQi 

Jo J4-4£^ Jo 



In the — limit, we have 



z^O Z 



Jo 



which hearkens back to the conclusion from Remark [(4)l to Lemma [LI] 

To confirm Eq. JST] simply differentiate both sides of Eq. [30] To prove Eq. [32] integrate Eq. [30] from to z; then 
Eg. [33] emerges as the specific case z = j for Eq. [32] To verify E as. [34] and [35] we combine the previously established 
results 

^ 1 /-iRCv^^) l + £ 1 i + sinyr 

G=— / log- — -d^=-/ log- ^dyr 

2jiJo i, 1-^ 4 Jo l-smi// 

with the cases z = ^ in and z - sini// Eq.[31] before carrying out trigonometric substitutions. ■ 

Remark (1) Our geometric derivation of the integral formulae in this proposition is partly inspired by Coxeter's in- 
vestigation of some functions relevant to non-Euclidean geometries 1.29.1 . in which context Catalan's constant G also 
appeared as the volume of an octahedron in hyperbolic space. 

An integral formula like Eq.[30]is a simple output of such a geometric approach. It is still possible to verify Ea.l30] 
combinatorially, using the moment expansion of elliptic integrals. This will be discussed in ^ 

(2) In [Re£ 6, Proposition 3.3], we used Ramanujan rotations to prove the following identity 

^^^'2 K(sin0)sin0d0 







1-X+ ^cos^e 



0<a:< 1. 



In Eq.[30] if we set t - cos^ 6 and z - s/l-x > for x £ [0, 1), then we obtain the identity 



iK(sin0)sin0d0 Li2(v'l-x)-Li2(-Vl-x) 



4 cos^ 



which bears close resemblance to the integral representation of [K(v^)]2. As a by-product, we have a strict inequality 
Li2(v^^)-Li2(-v^^)< v^^[K(v^)]2 for 0<x< 1. □ 

2.2 Dilogarithm Relations, Mehler-Dirichlet Decomposition, and Abel Transform 

The connection between the complete elliptic integral of the first kind K and the dilogarithm function Li2 (Proposi- 
tion l2.2D provides a wealth of integral identities of more or less geometric flavor. 



1 {l-z2)K(v^) zK(v^) , , 1 + 2 , , , 

df = — +logzlog-; , lz|<l, (36) 



Corollary 2.3 (Dilogarithm Relations) We have the integral identities 

rii^fS^d..r-j 

Jo 2\/4z2 + f(l-z2)2 Jo ^(1-22)2 + 422^ 4 1-Z 

^1 K(vT^) . riK(v^r^)sin0 ti^ 



r KVl-f r^K Vl-f sine TT^ / 0\ 71 

I dt+i\ — d^^ — +27:fl1ng ;:tan- , 0<e<-, 

Jo v/T^ Jo ^l-isin^e 2^2] 2 



(37) 
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-f 

Jo 



dt+ I — -c{t-2B]ctgtFm-, 
Jo 



\/si 



sm^ 9 - 1 



f*' K(v^) 



ji 

o<e<-, 

2 



/ C\t.-TT 

Jo V^^t 



arcsina;, 



0<x< 1, 



k(v/— ^|H^= TTsinVi ^x-Ti\os(x+ + 

VT+t\^^^t IV l + t) 



<X< +O0, 



the closed-form evaluations: 



L 
L 
I 



^ dt^---log^(^-l) 
2Vl+t o I 

log'' 

1 K(v^r^) , 71^ 1 



log2(v^-2), 

7^(3) 



2^4+7 12 6 

f [K(v^)]2df= f [K(v^)]2di „ , 
Jo Jo 2 

Jfi rjr2 T^fS)! r-*^ 

K(v^)K(v^l^)logidi= -71 — log2--^ = -tt / 
[ 2 4 J Jo 



4 



together with more integral representations of Catalan's constant G: 
7C(3) 1 r^KCr?)^ 1 



^r^^iogid, 

2?! ;7rJo 1 + 77 77 



G = -log2 + 
6 ^ 



2 3 

: + ^=^= + ■ 



G = -log(2 + 

O 



v'25-16f 4v/25^ ^625-576^ 
3 riK(v/r^) 



d^, 



/4-i 



-di, 



TT, 10+ V50-22v^ 5 



G = -log 



8 10- V50-22v^ 



5 />! 

- K(v^) 
o Jo 



\/G-2V^-t \/G + 2Vb-t 



At. 



(38) 

(39) 
(40) 

(41) 
(42) 

(43) 
(44) 
(45) 

(46) 
(47) 
(48) 
(49) 



Proof We note that the right-hand side of Eq. [30]is Li2(2;)-Li2(-2;) = 2x2{z), where Legendre's X2-function satisfies 
(see Eq. 1.67 in fs^) 
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1-Z 
T+2 



logz 1 + z 

+ I2(^)=y + — log — . 



This translates into Eq.[36l Setting z - jtan| in Eq. |36l we obtain Eq. [ST] In particular, when d - ji/2, both sides of 
Eg. 1371 vanish in an obvious manner. The imaginary part of Eq.[37l gives rise to Eq.[38l Meanwhile, as we read off the 
real part of Eq.[37j we arrive at Eq.[39l after the variable substitutions t^l-t, sin^O^ 1-x. (See also 4.21.1.2 in ifl^ 
for an equivalent form of Eq.[39l) By analytically continuing Eq.|39]to purely imaginary x, and noting that 

Jo Vl + tcos^e Jo y/l + t-tsin^e 1 + ^ 

we may confirm Eq. |40l 

The precise values of ^2(^-1), l2((^-l)/2) and X2(\/5-2) (see Eqs. 1.68, 1.69 and 1.70 in fsO*]) had been known 
to Landen, which correspond to the evaluations given in Eqs .l4lll42] and[43l respectively. To prove Eq.[44l we integrate 
Eq.[39]with the aid of the following formula (item 804.01 in Em) 



"'2 1 + sin 

log 



d0 



ib/a) l-sin0 Va2sin2 0-62 « 



= -K{b/a), a>b>0. 



This leads us to 



.fW)]^d..riogii^[r 

Jo Jo l-sm0 Jo 



K(v^) 



vsi 



:dt 



sm^0-t 



dB- 



nnl2 I 

n \ e\og- 
Jo 1 



l + sin0 _ 77rC(3) 
do - — - — , 



sinfl 
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where the last integral can be evaluated with Fourier series expansion. (An equivalent form of Eq.|44]has appeared in 
Remark 1 to Theorem 2 in f5j.) Before proving Ea. l45l we reformulate item 800.01 in as 



Jo 



arcsinv/f logsinfl KiyDlogt - nKiVT^) 

Hfl ■ , 



and proceed with the calculation 
''^K{Vi)logt-nK{VT^) 



Jo 



t - sin^ 



Jf7l/2 c 
Re 
Jo 



1 K(v^) 



dt 



t - sin^ 



Iogsin0d0 = 7rj [- - 0j logsinSdS 



a-[27r^log2 + 7C(3)] 
16 ■ 



Combined with Eqs.[T4]and|44l the formula above is equivalent to Eq.l45l 

The equivalence among the three quantities in Eq. |45] is not accidental, as the two extreme ends of Eq. |45] are 
actually connected by Ramanujan's rotation formulae [Ref 6, Eqs. 41 and 44]: 



[K(v^)]2= [ 
Js 



2i2-t)K(VX^ + Y^)da 
{2-t)^-fiX^ 4^' 



0<t<l; 



Js^ ^ ' 4:71 71 Jo 



1 + 1 



K 



d^, 



(50) 
(51) 



where {X,Y,Z) are Cartesian coordinates on the unit sphere : X^ + + - 1. Concretely speaking, we may use 
Eqs.[50]and[51]to convert Eg. [141 into 

^2 



2(2 - 1) 



S2 [ (2 - f)2 - fiX^ 



-1 



l-X^ da 



271 Jo 



K 



1 + k 



K 



1-k 



log- — —dk. 



This clearly supplies a geometric link between both ends of Ea.l45l 

To prove Eq.|46l we integrate both sides of Ea.[38lover the range 0<6< 7i/2. The left-hand side becomes 



JfTl/Z n 
Jo 

=-/'[/ 

Jo [J a 

= -[" f 

Jo [Jo 
Jo Jo J\ 





7l/2 



Sin 
d9 



' r K(V 
= dt+ I — — 

t Jo Vl 



1 K(v^^)sin0 



dt 



d9 



arcsin 
arcsin\/l-i 



KiVl^)dt + 



1 r^KiVT^), i + Vt 



2 Jo 



Vt 



log 



1-Vt 



dt 



de 



VT- 

du 



t - sin^ 



Jo 



K{Vl-t)dt+ Kh/l-r log 



1-^ 



d( 



y/l-il-t)u^Vl-u'^ 



K(\^T^)dt + 2 f ^^log-dT/=-f [K(\^T^)fdt + 2 f ^^log-dr] 
Jo 1 + V V Jo Jo 1 + n n 



after trigonometric substitutions and Landen's transformation. Meanwhile, the right-hand side evaluates to (see item 
1.442.2 in il2tl for the Fourier expansion) 



nil 12 

I 20 log tan 
Jo 



e r'^ 

-de = -4/ ( 
2 jo 



g cos(2^ + l)0 



e=o 



2£+l 



^ 2-(- 1)^(2^ + 1)71 7((3) 

d0 = 2 y ^ — „ = -27iG + 

^ (2^ + 1 3 2 



(52) 



f=0 



Quoting the result from Eq. |44l we thus have 



■iK(T?), 1 



-27rG + 7a3) = 2 f ^log-dr?, 

Jo 1 + ?? T] 

which implies Eq. |46l 

The inverse tangent integral Ti2 (z) :- -iX2(iz) has been treated in detail by Chapter 2 of llioll . where three functional 
relations 3Ti2(l)-2Ti2(5)-Ti2(^)- 5Ti2(f ) = f log2 [Ref 30, Eq. 2.28], 3Ti2(2 - n/3) = 2Ti2(l)+ f log(2 - v^) [ReflM 
Eq. 2.29] and Ti2(tan^)-Ti2(tan|f )+ |Ti2(l) = ^logtan^ - fglogtanfg [Ref 30, Eq. 2.51] paraphrase into our 
Eqs.|47l|48]and|Ml ■ 
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We may produce more integral formulae out of Eqs. [39] and |40l as explained separately in the following two corol- 
laries. 

Corollary 2.4 (Mehler-Dirichlet Decomposition) Let n £ {1,3,5,...} be an odd natural number, then we have the 
following integral formulae: 



j K|sin^|sin/3cos^d/3= < 0, 



2(n+l)/2(re-l)2 



(^)! 



n^l 

n = 1 mod 4 and n> 1 
n = 3 mod 4 



(53) 



which implies the identity 



1 K{k)rk\/l-k^ 7iK{\r\) 

■ak = 



1 + rjo 1- 



-2r(l-2^2) + ^2 ' 8 16(1 + r)' 
as well as two more integral representations for Catalan's constant G: 



-l<r<l 



G = I log 2 + ^ K ( v^l^j log V^l^d )f 
— f 



G = Zi(^-JLrK(sine) 



Li2(e^"')-Li2(e-2'S) 



I cost 



d0. 



(54) 

(55) 
(56) 



Proof We first set t - sir? ^,x- sin| in Eq.[39j and recall the Mehler-Dirichlet formula (Eq.[5ll: 



f 

Jo 



Y_ P/(cos0)cos 



(2£ + l)p 



K|sin-]sin/3d/3: 



Jo 



V2(cos/3-cos0) 2 



Meanwhile, the rightmost term in the equation above can also be developed into a series of Legendre polynomials (see 
8.925.1 of d) 



TiB n'^ \ ^ 4m -1 



■-l (2m - 1)2 



(2/n-l)!! 
2'"m! 



P2m-l(cOS0) 



As we compare the coefficients of P^(cos0), we may conclude that for positive integers m = 1,2, . . . , there are the follow- 
ing integral formulae 



JKlsin — I sin /3 cos 
I 2j 



(2^+l)/3 



= 

' = 2m 



71^ 4m-l [ (2m-l)!! l2 / - 2w - 1 
4 (2m-l)2 [ 2'"m! J ' * - J- 



which are equivalent to the statement in Eq. [53] 
Considering the familiar Fourier expansions 



~ ^ {2m- 1)15 Kl-r)cosf 

> r cos = IT 

^1 2 l-2rcos;3 + r2 



-l<r< 1. 



we may compute 



Jo 



K|sin^|sin/3 



r(l-r)cos ^ 
l-2rcos/3 + r 



■d/3 = 4 f^K(^)^vT^— 
Jo 1- 



r(l-r) 



2ra-2k^) + r^ 



dk 



71 

T 



^ 4m -1 
'■"^^i (2m -1)2 



(2m -1)!! 
2'"m! 



„2m 



7r(l-r2)K(|r|) ;r2(l-r) 



which entails Eq.| 
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As we integrate Eq. [54] over the range < r < 1 and exploit the fact that 2G = /o^K(r)dr, we may verify the equahty 



4:fo ^ 



K(^) arccos^ + 



4I^ld^ = ^-^log2. 
vT^I 4 16 



Next, we calculate (see also 4.21.4.8 in [18]) 



Jo 



K(^)arccos^d^ : 



Jo y/l". 



de 



Jo 



u^u , "^'21^12 (sine) -Li2(- sine) jt\ ^ 

arccosKd«= / ad- — log2, 

2sin0 4 



so that 



4 Jo 



. kiosk , , JiG TT^ , „ 

K{k)—=i= dk^ log2 

^/T^ 4 8 



rearranges into Eq. [55] after the variable substitution k ^ V1-k^. Recalling Eq.[46] we multiply both sides of Eq.[ 
by log J , then integrate over r e [0, 1], to deduce 



nnl2 



K(sin0) 



Li2(e^"^)-Li2(e-^"^) log2 
16jcos0 4 

^7I/2^ 



sin 6 



dO.'-l'-^ 
2n 



G 



n 

— log2. 
16 * 



This leads to Eq.[56]after we recall from Eq.[T]that /o*^K(sin0)sin0d0 = 
Remark It is not hard to recognize that we may deduce Eq. [53] from Tricomi's Fourier expansion 



K(|sin0|)= £ 

ra=0 



r(w+i) 



sin(4«+l)0, O<0<7r 



and vice versa. 

Corollary 2.5 (Abel Transform) (a) Define the Abel transform as 

2f(r)r 



(57) 

□ 



(Af)(. 



x):= [ 

Jx 



dr, 



then we have the following integral identity 



r 

Jo 



' x(Af){x) ( X 



K 



+ x^ 



dx-n \ f(r)rlog(r + Vl + r2)dr 
Jo 



(58) 



so long as the function fir) guarantees convergence of the integrals on both sides. Consequently, we have the integral 
formulae: 



Jo 1 



K(^) 



Vk^ + aHl-k^) 



arcsma 



-(TT-arcsina), 0<a<l 



dk^ < 



71^ log^(Va^ -1 + a) 

— + , a>l 

8 2 



a 

Jo 



kKik) 



(1-/^2)2 

1 kKik) 



Va^ - (1 + a^)k^ dk^^ \il + 2a^)logia + V 1 + a"^)- aV 1 + a^] , a > 0, 

o L J 



(1-/^2)3- 



71 

8 

[a'^-il + a'^)k'^f'^dk^ ^ [(8a^ + 8a^ + 3)log(a+ Vl + a2)-3a(l + 2a^)Vl + a2l , a > 0, (61) 
128 ^ J 



(59) 



(60) 



and 



kKik) 

(1-/^2)2 



r^— ^2 ^a^-il + a^)k^+aVT^ 

aw a^ -il + a'^)k'^ 1or 



[2 1 n 2rc r 2^ a(3 + 2a2) r T 

^n\-- —V 1 + a2(8 + 5a ) + — log(a + V 1 + a2) 

[ 9 36 12 



k 

, a>0. 



dk 



(62) 
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Furthermore, we have integral representations ofG and n /9>: 



G 



1 



_ 1 p fetanh 
TT Jo 1 



71^ 1 arccos/ 



''-K(k)dk, 



K{k)dk, 



along with two formulae valid for 0<6 < 7i/2: 



I ri ^tanh ^ Y J 



-f 

71 Jo 

if 

71 Jo 



a-k^)\/k^ + (l-k^)sec^0 



K{k)dk--\ logtan 

2 1 2 ® 4 



tanh ^Vl-k^ log cost 



(b) We have a pair of integral formulae: 



:tanh 



4 



which lead to the relations 



2[E(v^)-(l-u)K(v^)]= rK{y/i)dt^ /"' arcsinV? 0<u<l, 

Jo Jo \/u-t 

--[' ^Y.iJ^yt-j''^^'^'^'^'' 

Jo ^TT7 IV \ + t] Jo 



dt, ^>0, 



L 



'^E(Vt)-a-tmVt)^^ 7t^ n- ^ , ■ ^-^ a 

^ dt - — [vu(l-u)-(l-2u)arcsinv'«J, 0< u < 1, 

^Ju-t 4 

Proof (a) We first rewrite Eq.l40las 

7rlog(r + x/T+T^) = r 
Jo 



2a: 



-'0 \/l + x/r^ - V n/I + x^ 
then multiply both sides by /■(r)r and integrate over the positive real axis, as follows 

2x 



dx, 0<r<+oo, 



noQ noo r nr 

71 i f{r)rlog{r+Vl + r'^)dr- i f(r)r\i 
Jo Jo [Jo 



=K 



2f(r)r 



Jo V^l + x2 V Vl + X^i [a 



Vr^ -x^ 



dr 



dx dr 

r 

dx- I 

Jo 



'xiAf){x)^l X 



+ x^ 



(63) 
(64) 



(65) 
(66) 



(67) 
(68) 

(69) 
(70) 



dx, 



hence the formula in Eg. 1581 

Taking advantage of the table of Abel transforms (Table 13.9 in we may put down the following equalities 
involving a non-negative parameter a > 0: 



X /I 1 

Jo y/T+^\x 



K 



f 

Jo 



+ x^ ' V vTTx^ 

°- 2xVa^-x^ ^1 X 



+ x^ 



f 

Jo 



«4x(a2-x2)3/2^/ x 



3^1 + x2 Iv^ + x 



Jo v^^^ 



- x2 - x^ cosh ^ — I K 



x^ 



dx^7t[ \\ 5-]rlog(r+ \/l + r2)df 

Jo I H a'^ + r'^ I 

dx-7i i rlog(r+ \/ l + r^)dr, 
Jo 

dx-7i [ (a^-r^)rlog(r+ vTTr2)dr, 
Jo 

ra 

dx = TT / {a-r)rlog(r+\/l + r^)dr. 
Jo 



As we perform the variable substitution x ^ k/V on the left, and complete the integration on the right, we obtain 
the integral formulae listed in Egs. 159162] One might note that the a-1 specific case in Eq. [59] corresponds to Eq.[l] 
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We note that Catalan's constant can be cast into an integral in the following manner (of. 3.521.2 in ifli): 

'^dy 1 ydy 1 r°° log{r+ Vl + r^}dr 



„to(2Ai+l)2 „tb io Jo l + e-2y 2 Jo coshy 2 Jo 



1 + 



(71) 



where the last step involves a substitution y - log(r + Vl + r^). Meanwhile, the Abel transform reveals 



f 

Jx 



dt 



which brings us to 



f 

Jo 



oo xtanh""^ —ri: 



1 , 1 1 

tanh"^ 



^iT^ 



1 + 



K 



X ] C"^ X SI] 

s/l + ^l Jo 1 



o° X sinh ^ i 



K 



+ x^ 



+ x^ 



dx - nG. 



After replacing x with k/Vl-k^, we obtain Eq.[63]as claimed. (It is easy to recognize that Eq.[63]becomes Eq.|2] after 
the substitution k - Vl-k^ and Landen's transformation.) The proof of Eq.[64]is entirely similar, and is based on the 
facts that (cf 3.521.1 in IH) 



1 ydy 1 r°° \og(r + Vl + r^)dr 



n^^l r ydy ^ 1 r 
8 2 Jo sinhy 2 Jo 



rVl + r^ 



and 



1 

Jx rVT+ 



dt 



- arctan 



-• f 

X Jo 



°° arctan - 



K 



dx - n — . 



+ J.2 Vr2-x2 X X Jo Vl + x^ Wl + x^^ 

The foregoing approach naturally extends to the scenario where 

71-26] i cos0[Li2 iie-'^)- Lis (- ie"'^)] 



(72) 



(73) 



1 r°°log{r+Vl + r^)dr Bcoseiin 
2J0 sec^e + r^ -^[y-^^g^'^"- 



n 

O<0< -, 
2 



and 



Jx 



dt 



cos B sinh 



-1 1 



cos0tanh ^ 



\/jc^cos^0 + l 



}x sec^ + _ ^2 Tsec^w + x^ Vsec^y + x^ 

This leads to the formula in Eq.[65l (In contrast, there does not seem to be an "elementary" generalization of Eq.[64| if 
we have (a + rWl + r^ or + in place of + in the denominator, the resulting formula will involve the Meijer 
G-function, instead of the dilogarithm.) 



Multiplying both sides of Eq.|65]by sin0 before integrating in 6, we obtain 

-r 

Tl Jo 



tanh-1 v^l - ^2 logcose \/l-k'^s\n^B 
k{l-k^) kVl^ kii-m 



tanh"^ 



'K{k)dk 



-G+ — -cos6\ — I — - logtan 



n-2B\ i[Li2(ie-''*)-Li2(-je-''')] 



As we multiply both sides of Eq.[65]by cos0 and add to the equation above, we have 

-r 

n Jo 



tanh ^ Vl-k^ log 



cost 



:tanh"^ 



We will arrive at Eq. [66] after eliminating G from the last formula using Eq 



l-k^sin^B 



K{k)dk^G+—. 

4 



(b) The leftmost equalities in Eqs. [67] and [68] can be verified by differentiation in u and ^, respectively. Meanwhile, we 
note that the solution to the Abel integral equation (which is morally equivalent to the Abel transform) 



Jo 



^ cl>iy) 
Vx-y 



dy- fix) 



(74) 
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satisfies [Ref.lli, p. 39] 

Jo Jo Vx-y 

Specializing this to Eqs.[39]and|40l we obtain the rightmost equahties in Eqs.[67land[68l 

Clearly, the extreme ends of Eqs. [67] and [68] represent another pair of Abel integral equations. Their respective 
solutions bring us 

f''E{^/t)-{l-t)K{^/t) n r ^ i r- 

/ ^ at-—\ arcsmvidi = — [vu(l-u)-(l-2u)arcsm\/u], 0<u<l, 

Jo ^u-t 2 Jo 4 

as stated in Eqs. [69] and [70] ■ 

Remark The closed-form evaluations in Eqs. 159162] are not unexpected, as all of them can be regarded as appropriate 
integrations of some previously known identities. As one differentiates in the parameter a, one can reduce Eg. 1591 to 
Eq.[3ll Eq.[60]to Eq.[iO] while both Eq.[6l]and[62]become a certain multiple ofEa.[60l 
Differentiating Eq.[65]in 6 and multiplying by l/sin0, we obtain 



1 KiVT^) 



-( 

71 Jo 



cos^ 6 + sin^ f 



itanh"^- 



\/ cos^ 6 + siv? 6 V cos^ + k^ sin^ t 



dK=-^, o<0<5, m) 

sm20 2 



a result that will be derived again in SS.ll using a different method. Taking derivative of Eq. [66] with respect to 6 and 
dividing the result by sin20, we get the same relation as in Eq.[65]. □ 

3 Beltrami Rotations Revisited 

3.1 Beltrami Transformations and Some Simple Consequences 

In [Refl^ §3.1], we introduced the Beltrami rotations, which amount to the following integral identities; 



K(«=H'H^™1, 0<k<l. (75) 

ttJo 



Kir): 
K(T])-. 



n Jo l-f^(l-K^) 

2 fi (l + r)K(v^l-K2)d)f ^ 

- —z ,2.2 ■ 0<r<l, GSL) 

ttJo (l + rr-4r)c^ 

2 /■i(l-r?)K(v^^)dK „ -, fnms 

- —FT ^2 . 2 — ' 0<77<1. mi) 

n Jo {l-T]y + 4riK^ ^ 



As pointed out in (al, the four identities above are related to each other by standard transformations for complete 
elliptic integrals of the first kind: Eq. I75II descends from Eq. [75] and the imaginary modulus transformation; Eq. I75kl 
results from Eq. [75] and Landen's transformation; Eq. I75|^l is attributed to Eq. [75] along with a combination of the 
imaginary modulus and Landen's transformations. 

We have proved Eq. [75]in iQl using rotations on a unit sphere, which is inspired by Beltrami's work |32]. In the 
next two short paragraphs, we offer alternative derivations of the Beltrami transformations using the Mehler-Dirichlet 
decomposition and Abel transforms. 

For < r < 1, we may use the variable transformation k - VI-k^ to rewrite Eq. [54](a consequence of the Mehler- 
Dirichlet decomposition) as 

1 K(^l-K2)rjf2d)f 8 1 r^^.r / Z\ . 2 /■l(l + r)K(Vl-)c2)d)c 



K(.).^^ r K(/l^)rK2dK 8^ .1^ 2 r 

Til + rJo 1 - 2r(2)f2 - 1) + r2 7r4(l + r)jo ^ n Jo 



(l + r)2-4r)c2 

thereby demonstrating Eq. I75H (If we start with -1 < r < instead, then Eq.[54]leads to a verification of Eg. lTSp ) 
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As one combines the following identities 
' rlog(r + VTTr^)dr 



f 

Jo 



^Jl-em) and r 

Jx I 1 



rdr 



o<e< 1 



into the formula 



(76) 



one succeeds in proving Eg. lTSII using the Abel transform (Eq.[ 

At this moment, we may immediately appreciate the equivalence between two previously derived integral represen- 
tations of Catalan's constant 



G 



= - CK{k)dk^— f 
2 Jo 271 Jo 



iK(Vr^) 1 + K 
log- djc, 

K 1-K 



using the rational integral transformation given by Eq. (75] In other words, the Beltrami transformation (Eq. [75l l is a 
formal generalization of the "box-to-box transformation" (Eq.O introduced in Remark [(5)l of Lemma ll.ll Then, we may 
also deduce alternative forms of Eqs.[T]and[2]as follows: 



f'^Kik) riK(Vr^)/ 1 + k\ 
— - I d^ = / log hJflog IdK, 

8 Jo 1 + k Jo 7t{K^-l) [ ^ 4 ^l-Kl ' 



1 — i<-2) 

5^ [n^ - 3(1 + k)U2{k) - 3(1 - K)Li2(-)c)] dx. 



Gb) 
41b) 



Hereafter, an equation numbered xb is derived from Eq. x and a Beltrami transformation (Eq. [75l l. Descendants of 
Eq. X after transformations in Eqs. l75lll75Kl and l75£l will be labeled accordingly as xib, xlb and xilb- We will also refer 
to the four avatars of Beltrami transformations as types "B", "iB", "LB" and "iLB" for short. 
Applying the identity 



K(v^l-?c2) 



[^2cos2(/3/2)-Hsin2(/3/2)])c2 



to Eq. m we arrive at the expression below: 



1 rP 71 -a 1 K(vT^) _-, 

G+ — I logtan da=— / — tanh 

'^Jo 4 7iJo K\/l-K^sm^{B/2) 



KC0s(/3/2) 



■)c2sin2(/3/2) Vl-K^smHl3/2) 
After we differentiate both sides of the equation above in fi, the integral formula 

)ccos(/3/2) 



dx, 0</3<jr. 



et) 



1, 71-15 sin(/3/2) KiVl-x"^) 

- logtan / T 

4 4 271 Jo l-K2sin2(/3/2) 



Kcos(i3/2) ^ , 
:tanh 



Vl - Jc2 sin2(/3/2) Vl - Jc2 sin2(/3/2) 



d)c, 0<jS<7r igL) 



arises as a result0 Here, in Eq.^J the bracketed term is an even function of jccos ^(1 - sin^ ^)~^''^, so its value is 
not affected by the choice of the univalent branches of the square root. By anal3rtic continuation, Eq.^] extends to the 
following identity valid in the open unit disk: 



1 1-z _ z(l + z2) K(V1-k'^) 

4 + ^ ~ n Jo (1-Hz2)2_4^2^2 



7f(l-z2) 



V^d-H 22)2 _ 422^2 

If we set z - itan(0/2),O < 9 < 7i/2 in the equation above, we recover Eq 



tanh"^ 



v/(l + z2)2_ 422^2 



dx, \z\<l. 



^Of course, one can rewrite the right-hand side of Eg.B^bv the substitution 

[ 21 nJo 1 



2 fl K(Vl-K^)dK 
k2 sin2(/g/2)' 



0£P<JI. 



A similar procedure can be performed on a few formulae in the rest of 33. li as well, except that some subtle changes in the domain of validity may 
occur. The details about paraphrasing these integral identities will be omitted in 33.11 but will appear in the collation of formulae in Appendix IaI 
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In Remark [(2)1 to Proposition 11.21 we mentioned Eg. [TTl without a proof. Now, we can demonstrate Eq.[l7]by con- 
necting it to Eqs. [Hand [55] via the type-B transformation (Ea.[75ll: 



/;H-fiT^fii 



1 2K(vT^) / 1 



1 djc 



— = f^KfVl-)f2)log\/l-K2d)f =-2G + 7rlog2. 

k n Jo y I 



Setting < < 1 in Eq.[30j writing k for \/l-t, and applying the rational integral transform in Eq.[75l we obtain 
an equivalent form: 



1 4K(vT^) 



>'0 7rK\/(l 



1 , V^(1 + z2)2k2-4z2 + (1-22)^ 

log —= + log ^ 

vl-?c2 + z'^)'^k'^ -Az"^ + {1 + z'^)k 



6-K ■■ 



Li2 (0)-Li2(-0) 



dsoh) 



which extends to be valid for all the points satisfying < 1, by analjrtic continuation and the Vitali convergence 
theorem. In the z ^ limit, we obtain the identity 



ri K(^l^) ^^^ 1 

Jo 1 - ^2 



A.K - n. 



(77) 



Differentiating Eg. lSOhl in z, we arrive at an equivalent form of Eq.[3T|for \z\ < l,lmz ^ 0: 



Jo 7t 



4k(1-z4)K(vT^) 



[(1 + 02)2^2 _ 4z2]^(l + z2)2^2_422 



1 Vil + Z^)^K^-4:Z^ + (1-Z^)K 
log —= + log ^ 

v1-k:2 + z^)^k^ -4:Z^ + (1 + z^)k 



djc 



1 802K(v^l-)f2) 



-r- 

Jo yi\ 



1, 1 + 



■dK - -log- 



r[(l + 02)2^2_422] ■■ 2"° 1-2' 

Similarly, Eq.[32]has a cousin 



([3Tb) 



L 



iK(Vl-)f2) 
47rx:2 



log2(l-)c2)-4 



1 V(l + z2)2^2 -Az'^ + {l-z^)K 

log , + log , 

V^l-)f2 + 22)2)^2- 4^2 + (l + z2)^ 



dK 



= 2Li2 (Z) - 2Li2 (-2) - (1 + Z)l0g(l + 2) - (1 - Z)l0g(l - Z), |Z| < 1. 

In particular, the special case oiz-i for the equation above leads to an analog of Eq. 



1 riK(Vl^) 

log:^ - - 
4 2^ 271 Jo K' 
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logd - jf)log(l + jf)djf. 



([32b) 



([33b) 



The following identities also follow naturally from the Beltrami transformation: 
4 f ^ K( v^l - k2 ) arcsin( jf x) 



4 r 

\x - — \ 

71 Jo 



1 . A 2^ 4 rK(V 

jrlog(x+ v l + x^)= - I 

71 Jo 1 



4 r'^KiVT^) 



■dx. 



arcsin 



: arcsin 



Vl + X^l Vl + (l-)f2)x2 Wl + x2 



x£[0,l], ([39b) 
djc, a;e[0,+oo). ( |40b ) 



(Here, an equivalent form of Eq. [390 appeared as 4.21.4.5 in 118].) If we analytically continue Eq. |39fa| from real-valued 
X to purely imaginary x, we obtain yet another integral formula (see also 4.21.3.10 in lEi): 



7rlog(a; + Vl + x^) -- [ 
71 Jo 



4 riK(v^l-K2)log(Kx+ 



dx, 



kVI + k^x^ 

When a similar service is performed on Eq. l40fal we may deduce 



xe[0,+oo). mt) 



4 riK(Vl^) 

7rarcsina;= — / ^ — 

7r Jo 1 - K"^ 



sinh ^ 



Vl-x'^l ^1-(1-Jc2)x2 



:sinh ^ 



KX 



dJc, x£[0,l) ([40E) 



instead. 
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With the Beltrami transformation in Eq. [75l one may readily verify the following integrals ofK{Vl-K^) times an 
elementary function of k, which may otherwise appear formidable: 



L 



1 2KiVl - k^)[k(1 + r)v^tanh-^ sfF- rVl - arcsinjc] 7rK(|r|) tt^ , , 

djf= , 054h ) 

;r(l + r))f[4r + (l-r)2jc2] 8 16(1 + r) 



2 r^KiVl-K^) 
n Jo 



log(l + a)-)ctanh ^k- 



\/l-(l-jc2)a2 



log 



Vl-xHaK + v^l-(l-K2)a2] 



■(TT-arcsina), 0<a<l 



7r2 log2(v^^2Tl + a) 
r + :: , a>l 



and 



2 r'^KiVT^) 



2 r^K(v 

71 Jo 1 



(l + a2 + K:2-a2)c2)arctana a K\/l + il-K^)a^ 



2(1 -jc2) 2 
= I [(l + 2a2)log(a+ \/l + a^)-a\/l + a'^j , a > 0. 



l-jc2 



arctan 



aK 



V'l + (l-)c2)a2 



djf 



((Bob) 



The output of the type-B transformation (Eg.lTSll always takes on the shape of /q^K(\/1 - K^)pe(K)AK where pe(f) is 
an even analytic function of jc. However, the type-iB transformation (Ea. l75lt 



dxr 



results in a host of integral formulae for /q^K(v^1 - k'^)p(k)Ak in a different flavor, where p()c) is not necessarily an even 
function in k. For example, we can transform a formerly derived integral representation of Catalan's constant into 



2jo 2jo 1 + Jf 2jo l+x/w2x/w2 ;rjo ^ ^ 



KarccosK 



v^l^ 



■log(2)c) 



d)c, (78) 



after two rounds of applications of the type-iB transformation (Eq. |75l| l. 

We gather below some other consequences of the integral transform in Eq. |75ll while leaving out the computational 
details of their straightforward proofs: 



n 
G 



riK(0 1 r'^ - 

---^d^=- K(sin0) de 

Jo l + i 71 Jo cosy 

71- a ^ 1 r'^KiVl-K^) ' 



4 Jo 4 
2G 1 r'^KiVl-K^) ' 



,71 sine -26 

COS0 

logtan da + - I 5 

"Jo 1 - Jc^ 



2G _ 1 r K(V 
T~2jo ~ 



2. 

2k 



Vcos2(/3/2) + K2sin2(/3/2) 



dx, 



2G_1 r^K(Vl^) 
~5~~2jo ~ 



Vi + 3k2 

(/S + Dk 

Vl + (5 + 2^5)^2 Vl + (5-2v^))c2 



(v^-l)K 



7^(3) 1 K(v^l - jf2)[7r - 2kK(v^1 - k^)] 



'Jo 



1-Jf2 



dK, 



-logtan 
4 4 



Li2(z)-Li2(-z): 



p K(Vl- 
Jo ^^2 

7r-/3 sin;6 



ji 
2 



^ KiVl-K^) KarccosK 



dx, 



-! 

Jo 



kK(V1-k2) 



8 Jo [cos2(/3/2) + )c2sin2(/3/2)]3/2 
2 r'^KiVl-K^) 



dx, 0<p<7i, 



GIb) 

djc, 0<p<7i, (|1b) 
(|lB,with;S = 2jr/3) 
d K, (Hb , with p = 2jr/5, 47r/5) 

iSIb, integrated over < /3 < jr) 
(IHb X cos ^, integrated over < /3 < jt) 
(IHb, derivative in jS) 



2 r^K(Vl-i 
jr Jo 1 - jf 2 



1 

log— — + 

1 + 2; 
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2ZK 



40(1 -z2) 



log 



n 

T ' 

1 + z 



V(1 + z2)2k2-(1-z2)2 

2 riK(vT^), 

-/ ; 2 logJfdK, 

Jo 1 - 



: tanh 



.1 V(1 + z2)2k2- (1-22)2 



(1 + z2)k 



dK, |z|<landz^7^1, ((318) 
(|30]r. withz^: 1) 



4 K(\/1 - k2) / TcarccosK 
Jo I - »^2 



: p K(Vl- 
• Jo 1 - f 

/' 

Jo 



(l + z2))f 



^(1 + z2)2^2_(i_22)2 



d)f. 



: tanh 



.1 v/(l + z2)2y2_(i_g2)2 
(1 + z2)k 



IIb, derivative at z = 0) 
K(Vl-)c2)d)c 



(1 + z2)2k2-(1-z2)2' 



IzKl, 

(HIIb) 



Trarcsmx : 



if 

TI Jo 

4 

jr Jo 



tanh X - 



V'1-(1-k2)x2 



: tanh 



1-- 



V1-(1-)c2)x2 



tanh"^ 



Vl-(l-)c2)x2 



K(vT^)d)c 



l-jc2 



, 0<a;<l, 



V^1-(1-)c2)x2 



K(Vl-Jc2)dJc 



1-(1-k2)x' 



2W2 ' 



0<x< 1, 



;rlog(a;+ \/lTx2)= — f tanh ^ 
jr Jo 

Jl Jo 



jcx K(\/1- Jf2)djc 



< X < +00, 

KX . , _i KX \ K.(V 1 - K2)d)C 



\/1 + k2x2 )C\/ 1 + )f2x2 



tanh"^ ■ 



V1 + k'^x^I 1 + )c2x2 



< X < +00. 



391b) 



B, derivative in x) 
(ESIb) 

(|40]r. derivative in x) 



As we make the substitution x '-^ x/VTTx2 in the antepenultimate formula, we also obtain the last formula. It is also 
worth noting that Eas.l39 ^l and l40lTtl are in fact identical. The same can be said for the pair of equations 140 ^1 and I39Ir I 
Applying the type-LB and type-iLB transformations to Eq.[l] we obtain the following identities 



8 Jo 1 + r Jo 

^.f!^d,.r 

8 jo 1 + r? Jo 



1 2K(V1^) 



Jo 



dr 



(l + r)2-47-Jc2 



Jo 



iK(va^) 



arcsinjcdjc, 



i2K(v^l-?c2) r /■! 1 (l-ry)dr? 



11 1 



/■iK(v^l^), ^ 

dK=-/ — ^logJcdK 

Jo 7r(l-K^) 



+ (1 - r;)2 + 4riK^ 



•EIlb) 



Neither of them is particularly new, as Eg. llt Ri follows from the x = 1 specific case in Ea. l39fal or a variable transforma- 
tion k - Vl-K^ in Eq.[64l and the z ^ 1~ limit of Ea. [30lR] recovers Eq. IIItrI Doubtlessly, the two equations above also 
arise from the applications of type-B and type-iB transformations to the form of Eq.[T] after Landen's transformation: 



8 2 Jo I ''Jo 7t [Jo l-(l-^2)^2 

'iK(Vi-)c2) r 



d^ 



1k= f 
Jo 



71K 



VT- 



arcsin)cdK, 



^arK(vCT)d^.r^^^Mr idK^-f^^^^og^dK. 

8 2jo I -' Jo [Jo l-(l-.^2)(i_^2)J ;j(l_^2) 

In the next subsection, we will perform the transformations specified by Eos. l75land l75ll on multiple elliptic integrals 
which have undergone Landen's transformation beforehand. In these contexts, we will be effectively using the trans- 
formations of types "LB" and "iLB" in the disguise of types "B" and "iB". 

3.2 Duality Relations 

The next corollary presents some further applications of the Beltrami transformations (re)introduced in 33.11 
Corollary 3.1 (Duality Relations) (a) We have the following identities for r > 0: 



W(r) :-- 



'irK(Ode riK(vT^)d)c 
v^lT72 



/•I rK(Od( _ 
Jo r2+^2 I 

2 

n Jo 



. + r^ + Kr 



KrarccosK / r 

— _ V 1 + H logjc 

Vl- " 



J 1 + /-2-K2 ^1772 



K 



(79) 
(80) 
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2 I Kx/l + r^arcsiiiK / r 

— I , hrlogv l-K^ 



= -j^ rlog 3—-^ + —^^^=^K\-^==\. (82) 



This implies two more integral representations of Catalan's constant 
as u;eZZ as the functional equation for Wir),r> 



(83) 



W(r) + j W ( j V^l + r2] = 4j (r + v^l + r2] W |i (r + v^l + r2]^| (84) 
rirK(£)d£ ri vT+72K(ade _ 4{r+VT+7^fK{T])dTi 



(b) We have the identities 



M,r, :- /; = -i jVlog VIT^^H^™ . J£I^K| ^ | (85, 



Vl + r^ \ Vl + r^ 
K(Vl-K2)d?c log(r+v^l + r2) f 1 \ 
JiJo 1 + r^-K^ lv^l+72J' 

^/le integral formula 



/."( 



r2 + )f 2 1 + r2 - )f 2 



71 log: 



Kfx/rr^WrHr^-— — !— (87) 



a/id ^/le functional relation 



W(r) + M(r)= ^ — 7r+\ 57^ = - ''log 9 5 = - Vl + r2log- 

Jo H + f ^ Jo 1 + Jo K H + JC^ TT Jo 



1 + r2 - k2 



" ' Kf^l (88) 



2 \/l + r2 V \/l + r2 / 
all valid for r>Q. As a consequence, the following integral relations hold for x > 0; 



_i /l + x2jc2K(\/l-?c2)dK fi _i /x2 + ?c2 K(Vr^)dJc r 

I tanh T + / tanh M/ ^ , — - I 

Jo V 1 + ^^ Vi + x2k2 Jo V 1 + ^^ Vx2 + k2 Jo 



1 /l + x2jc2 K(\/l-?c2)dK fl , _i / x2 + ?c2 K(v^l-K2)dJc _i a; K(^l^)djc 
*' ' ' * — ^ *' - ' tanh ^ — , (89) 

\/x2+lc2 \/x2 + )c2 



K(v'l-?c2)d)c , -1 a: K(v^l-)c2)d)c 



Jo \/1 + x2jc2 \/l + x2)f2 Jo 



\/x2 + Jf 2 Vx^ + K^ 

and there are several multiple elliptic integral representations for K(\/T)K(\/1 - A),0 < A < 1; 

J(i-A)/(i+A) vT^T^ J(i + A)2jfe2 - (1 - A)2 Jo v^4A + (l-A)2^2 Jo 



(90) 



i(l-A)/(l+A) vT^T^ ^(1 + A)2jfe2 - (1 - A)2 Jo v^4A + (l-A)2^2 Jo v^(l-A)2 + 4A^2 

2 _i 2v^ K(v^l-)f2)d)c 

= — / tanh ■ 
jr Jo 



2 f"*^ 

= - / tanh"^ —— 
jr Jo ^ 



V4A + (1-A)2jc2 V4A + (1-A)2k2 
1-A K(vT^)d)c 



- A)2 + 4Ak2 V(1 - A)2 + 4Ak2 
Proof (a) With the type-iB transformation (Eq. I75II I and the elementary integration 



(91) 



2 r VT^d( 1 

r > 



TT Jo 



'0 r2 + ^2i_^2(i_^2) y/i77:2 + ^^' 
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one may reveal the truth of Eq.[79l Applying a second round of type-iB transformation to the right-hand side of Ea. [79l 
we obtain 



r^K(Vl-K^)dK _ ri K(0 fd^ _ 2 jcrarccosTf ^/j-^j^^K(r+ >/l + r2) 

Jo Vl + r^ + Kr Jo Vl + r2+rx/l-£2 n Jo V1-k^ Vl + r^ 



As we recognize that 



K(Vr^)djf 

1 + -K^ 



2 f^K{Vl^)dK 1 



2 pK(V 

71 Jo 1- 



- jf 2 1 + r2 1, vTTr 
we may confirm Eq.[80l Now, writing | - arcsinjc for arccosK, and noticing that 



ri Kr K(vT^)d)c _ rK(Od^ 
io ^fr^f l + r2-K2 "Jo r2 + e2 



(92) 



follows from a simple substitution k ^ we see that Eqs.[79]and[80]entail the identity 

/■i Krarcsin)cK(v^l-)c2)d)f _ , ^ 1 K(Vl-Jf2)dJc Jrl°g(l+7^) / 1 

io l + r2-K2 "Jo l + r2-K2 2 vTTT^ 

By analytic continuation in r, we may convert the last identity into 

kVI + r2 arcsiuK K( Vl - jc2)d)c , lK(v^l^)d)c tt log(l+ "^^^j^/ 1 \ „^ 

Jo ^Tr^^ ^^^-\o ''""K^^- 2 ^ ^ItiTT^)- 

On the other hand, if we apply type-B transformation (Eq. [75l > to the right-hand side of Eq. [79l we may set up the 
following equation 



iK(v^l-?c2)d)c K(^) kdk 



riK(v^l-?c2)d)c _ 
Jo \/l + 7-2 + )f r Jo 



\/l + 7-2 + Kr Jo V^l + r2 +)C\/l-^2 x/l_^2 



2 

n Jo 



K Vl + r^ ar csin jc v^l - k^( r + Vl + r^) 

+ riog 



Vl-K^ Vl + r 

which implies Eq. [81] after we exploit the identity 



KiVl-K^)dK 



2 r'^rK(Vl-K^)dK 1 

: xV 



2 r 

Jo 



vTTr2 I vTT72 

As we set Eg. 1921 1 into Eq.|81] we obtain Eq.[82l In passage to the limit r +oo, Eqs. [79l and [Ml recover the identities 
stated in Eq.[78l while Eqs.[81|and[8llead to Eq.[83l 

After applying Landen's transformation on the right-hand side of Eq. [79l we get 

2K(r?)dr? 



ri^;K(Od£ r\ 
Jo r2 + ^2 j„ 



v^l + r2 + r + ( Vl + r2 - r);/ 
Setting r = sinhs in the equation above, we obtain 

e-W(sinhs).e- f ^ ^^nhsK(Odf 2K(r^^ , , 0. 



sinh^s + ^2 Jo 1 + e 2s^ 
By analytic continuation from s to s + ^ . we have 

^e-Wacoshs).e-r'-°-^f^^y-f'^- - 0- 
Jo coshes -^2 Jo 1-e ^"77 

Adding up, we arrive at the following functional relation 



STX.. • u ^ ■ J SinhsKCOd^ C0ShsK(0dn 4K{T])dri ,.2sTX..■2.^ 

e''W(smhs) + fe''W(fcoshs) = e'' I o + I o ttI^I -r-'^ ^ 4:ie'''W(ie'''), 

[Jo sinh"^s + f2 Jo cosh"^s-^2 J Jo 1-e ^^r/^ 

which entails Eqs.[84|and[84| after the back substitution s = sinh~^r. 
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(b) By the transformations of types B (Eg.lTSll and iB (Ea. l75ll l, we may establish the integral representations 

'iKK(Od^ 2 /l-K2K(^l^)d)f 
-- / rlogW — 
JiJo V 1- 



Mir) :-- 



Jo 1 



2^2 



2 /jcrarccosK / r . , \ K(\/l - Jf2)d?c 

= I — , vl + Hsmh r — , 

which simplify into Eqs.|85]and[86l As we combine Eqs.[82]and[85j we obtain 



W(r) + M{r). 
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-f 
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(93) 



whereas Eqs. [80] and [86] together entail the relation 

71 Jo 

— / Vl + r^ log 
jr Jo 



W(r) + M(r): 



, lK(va^)d)c iogvT+72^^ 

1 + H log ^ z 1- K 



K 1 + /-2-K2 v^lT72 l^lT72 

vT+72 K(^ 



1 + r2 - k:2 
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(94) 



Equating the two expressions for W(r) + M(r) in Eqs.[93land[94l we have verified Eq.|87l 

To establish the last equality in Eq. [SS] we set r - k/Vl-k^ for < ^ < 1, and make use of the rightmost term in 
Eq.[79]to put down 



which adds up with 
Mir) = M 



IvT^^/ Jo i + kK 

i7-CK(0d^ K=^/T^ rir)cK(Vl-jf2)d)c 



- f !l 
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p rxH 
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L + r2^2 jg l + r2(l-)c2) 

into a form anticipated from the Beltrami transformation in Eg. 1751 

k 



I ^KK(v^l-Jf2)d?C 

Jo r^F^^2 



W(r) + M(r) 



= K(ft) = K 



As we may recall that the last term in Eq.[79]itself is a consequence of the Beltrami transformation of type iB (Eq.lTSi 
it might be appropriate to name the identity expressed by the two extreme ends of Eq. [88] 



i7-K(0d^ /■irfK(Od^ 71 1 



Jo r2 + ^2 1 



io r2 + ^2 1 + ^2^2 2^l772"^^i^ 

as the sesqui-Beltrami transformation. 

We perform Abel transform on both sides of the last equation: 

1 y;:2T3^[jo r2+^2 1 + ^2^2 J 2I Vr2-x2 lio 



df 



which results in 
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Jo Vx2 + )c2 Vx2 + jc2 

On the other hand, using the Beltrami transformation in Eg. 1751 we may verify 

df ^ 

2 Jo ,/^^ 



(95) 



■"^io ^/l + x2f2 Jo 



tanh 



1 /l + x2K2K(v^^^)d?c 



1 + x^ 



\/l + X' 



^"^io 



tanh 



1 /x2 + )c2 K(v^l^)d?c 
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so Eq.[89]is clearly true. As we trade x for llx in Eq.[89l we obtain 



Jo 



tanh 



i+x2 



i-?f^)d)c r 



tanh 



1 + ^2 vT+A2 



/■^ 1 

= X / tanh"^ —^=^= 

Jo Vl + X^K 
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which gives rise to Eq. |90l 

In [Ref.lii Proposition 3.3(b)], it was shown that 



K(\/A)K(vT^)= 

J(l-A)/(1 + 



K{k)dk 



holds for < A < 1/2, and the right-hand side of the equation above is left unchanged after a reflection A 1 - A. As we 
carry out the following computations for 6 > 1: 
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before invoking Eq. |95l we can prove that 
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which verifies Eq. [911 



4 Moment Expansion and Hypergeometric Summations 



Some of the integral identities developed in § ^1131 can be recast into infinite series, where each summand is ex- 
pressed combinatorially. This leads to some hypergeometric summations possibly of independent interest. 
For example, using the moments ofKlVl-x^) l33tl 



j\''K[Vl-K^]dK^^ 



r(^) 
r(^) 



we may reformulate some previously derived results as 
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One may double-check these series representations by evaluating them directly with the Wilf-Zeilberger hypergeo- 
metric summation, as implemented in commercially available mathematical software packages, such as Mathematica. 
Combining our series for 2nG - 7C(3) with the following result of Gosper js^l : 



2nG-'-^^t ' 



n=0 
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we may deduce two more sums: 
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Here, the two series representations of Catalan's constant 



4^0 2^1+1 
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3n + 2 



7r„to('i + l)(2ra+l) 



{2n+m 



provide no numerical advantage over the expression G - T.'^^Q(-l)"'{2n+ 1)"^, as the summands all have 0(l/n^) decay 
rate. 

We could improve the convergence rate by turning other integral representations of G into series. It is worth noting 
that Eq.|48] translates into Ramanujan's series representation for Catalan's constant G [14]: 



G = - log(2 + 



3rlK(^T^),. TT, /TT. 3 ^ (5)" 
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while Eq.|49] corresponds to Bradley's series acceleration formula 
TT, 10+ v^50-22v^ 5 
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Both the series of Ramanujan and Bradley converge exponentially fast, and thus have provided efficient algorithms to 
evaluate Catalan's constant numerically. 

Making wise use of combinatorial techniques, we may also verify Eq. |31] in its most general form. Suppose that 
1 1 - 0^1 > 2\z\, we may use Taylor expansion to compute 



Jo 



2(1 -z^) ~ / 2z (-2Tn\ 

[(l-z2)2 + 422^]^(l_ ^2)2 + 4^17 ^~ {l-Z^ft^yi-Z-^] (271+1)!! 



In the meantime, one can verify that 
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both satisfy the differential equation 
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which extends to the open unit disk |z| < 1 by anal3rtic continuation. One may compare this mechanical derivation to 
the geometric motivations that led to Eq. [301 

Using the moment expansion, we may express the function M{r) (Eq. |85| l in terms of generalized hypergeometric 
series: 



M(r): 
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which is also a special case of 2.16.5.1 in 112511 . By the relation 



W(r) + Mir) = - K 
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we have 



W(r). 
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which may be compared to a particular case of 2.16.5.10 in |25]. As we analytically continue the functional equation for 
W(r) (Ea.[84t to r = - j sin0,O <6 < nil, we obtain a hypergeometric identity: 
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Here, in the last step, we have used Landen's transformation and the inverse modulus transformation to deduce 
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l + jtanej 2COS0 2 



Setting - nlA in Eq. |84H and employing the special value 



l\^2 



we obtain the relation 
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Thanks to the Chowla-Selberg theory |35, 3^ 37, S^l which leads to the special values 
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I 121 I 12i 



12^ v^"l 12 i 2™ n 
the case of = :?r/12 for Eq. l84b-| also results in an exact evaluation 
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Similarly, setting 6 - 5nll2 would bring us 
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It would be unfair to go without mentioning the following summation formula in Ramanujan's notebooks [Ref.|3£ 
p. 153]: 
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nto (271 +1)2 cosh ^^^^±i^ ~ 4K(7M 
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where the prescribed domain of validity ensures that 2;K(\/A(z)) - iK{\/l- X{z)), with Kz) being the elliptic modular 
lambda function. Using the imaginary modulus transformation for complete elliptic integrals, one can verify that X{z) 
satisfies the following transformation laws 
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in the upper half-plane where Imz > 0. Meanwhile, Landen's transformation for complete elliptic integrals entails the 
following duplication formula for the elliptic modular lambda function (cf [Ref 40, §135]): 
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With the correspondence between trigonometric factors and elliptic modular lambda functions (cf Eas.l96landl97t 
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where zli>Q, 



and the related modular transformations for complete elliptic integrals of the first kind 

is/Jiz) 1 
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we can rewrite Ea. l84H using Ramanujan's summation formula: 
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The special cases of Eq. [84F] with 6 - 7i/4,7i/12,5n/12 can be converted into some h3rperbolic series of Ramanujan type: 
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although these formulae did not seem to show up in Ramanujan's works. By anal3^ic continuation and complex conju- 
gation, one can further generalize Ea. l84Kl into 
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The left-hand side of Eq. I84F.^ b lso has well-defined limits when z tends to certain rational numbers. For example, the 
case where z ^ 1 recovers the familiar series 2:°°_„ l/(2re + 1)2 = 7r2/8. 

The moment expansions involving the product of two elliptic integrals also paraphrase some of our integral formulae 
into hypergeometric identities. For instance, we may rewrite Eg. 1 141 as 
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after invoking the expression for /q^[K(\/1 - 7c^)]^jc^""'"-'^ djc from Eq. 18 in iQl- Similarly, according to the formula for 
foKix)K{VT^)x^''^^dx from Eq. 18 in @1, our Eq.|45]can be rephrased as 



^2n + l 



I 1 1 

2' 2' 



n + l.n + 1 
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Interestingly, every summand in the last line is actually a rational multiple of tt^ [cf. Ref.@, Lemma 1]. This statement 
can be proved by using the third-order ordinary differential equatioi|^ 



t{l-t) 



dfit) 
dt 



satisfied by fit) - K(v^)K(v^l-i), which leads to the recurrence relations 
K(a;)K( v^l - x2] x^"'^'^dx = K{\/t)K{VT^)t'' dt 
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= 2 f^K(v^)K(^r^)f"-'idf-2 Cf'-^ltil-t}-^ 
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Jo 
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after integration by parts. Here, we have /q^ K( \/F)K( V 1 - i ) d i - n^lS according to Ea.[26land the recursion relations 
outlined above reveal the rationality of 7r~^/o'^K(v^)K(\/l-f)i"df for all positive integers n. 



^As one may recognize, this third-order differential equation just paraphrases [Ref 6, Eq. 16] in the case where v = —1/2. The use of such a 
Legendre-Appell differential equation dispenses the combinatorial argument in the proof of Lemma 1 in Q). 



40 



A Index of Multiple Elliptic Integrals and Related Series 



A.1 Brief Lists of Identities for G, ((3) and n 

We compile here all the formulae leading to Catalan's constant G, Apery's constant ((3), and the circle constant 
71, which had been derived somewhere in the main text. We keep track of the original equation numbers, followed by 
contextual annotations if necessary. Some closely related representations are (re)grouped together, so the enumeration 
does not strictly follow the original order of appearance. 

The versatility of integral and series representations for Catalan's constant has been well recognized, as early as 
in Catalan's own treatises 10, B- As we have biased towards results that have geometric interpretations and/or have 
connections to elliptic integrals, the list below is by no means an exhaustive catalog of identities for Catalan's constant 
G. In fact, we will forgo many definite integrals of elementary functions giving rise to G, which can be looked up in 
well-indexed references such as fl2^. Furthermore, a few formulae gathered below have been known long before the 
current work, and our "geometric proofs" for them do not necessarily represent the most straightforward approach. 
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Continuing the list above, we bring together some integral and series representations of Apery's constant ^(3), which 
do not explicitly involve Catalan's constant G. 
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The last part of this subsection contains a collection of formulae that evaluate to tt^/S, up to an integer power of Ji 
times a rational number. 
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A.2 Miscellaneous Integrals and Series Concerning Elliptic Integrals 

In this subsection, the results will be divided into two categories: constants and functions. The arrangement of 
equations within each category roughly follows their original order of appearance. 

We first present a set of identities that invoke special values of surds, logarithms, elliptic integrals, and (generalized) 
hypergeometric series. 
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Lastly, we provide a collation of integrals and series that evaluate to elementary and special functions. 
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